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I.  INTRODUCTION 


Research  on  microwave  generation  from  rotating  electron  beams  in 
various  conducting  boundary  systems  has  been  pursued  at  the  University  of 
Maryland  under  AFOSR  sponsorship  for  several  years.  Recently,  the  major 
emphasis  of  this  research  has  been  on  the  production  of  radiation  from 
rotating  electron  beams  in  magnetron-like  conducting  boundary  systems. 

These  experimental  configurations  are  now  recognized  as  an  entirely  new 
type  of  microwave  tube,  referred  to  in  the  literature  as  a  Gyromagnetron , 
High  Harmonic  Gyrotron,  or  Cusp  Injected  Magnetron  (Cusptron).  The 
interest  in  this  new  device  has  centered  around  its  potential  to  reduce  the 
required  magnetic  field  in  microwave  tubes  by  an  older  of  magnitude  by 
allowing  operation  at  a  high  harmonic  of  the  electron  cyclotron  frequency. 

This  area  of  research  has  seen  substantial  progress  during  the  last 
year,  and  a  summary  of  theoretical  and  experimental  research  conducted 
during  this  period  is  presented  in  Section  II,  Results  from  research  on 
millimeter  and  submillimeter  waves  produced  by  rotating  electron  beams  in 
rippled  magnetic  fields,  partially  funded  by  AFOSR  (with  additional  funds 
from  DOE)  are  also  summarized  in  Section  II  of  this  progress  report. 


II.  FINAL  PROGRESS  REPORT 


For  the  Period  December  1,  1984  to  November  30,  1985 

A.  Experimental  Research 

1 .  Improvements  to  the  Experimental  Facilities 

During  the  past  year  the  experimental  facilities  have  been  enhanced  by 
the  installation  of  a  completely  digital  fast  data  acquisition  system 
funded  by  a  grant  from  AFOSR  under  the  DOD-University  Instrumentation 
Program.  This  new  system  includes  five  channels  of  Tektronix  7912AD  fast 
digitization  controlled  by  a  DEC  11/73  computer.  Additional  equipment 
related  to  this  project  includes  an  Imagen  laser  printer  and  several 
graphics  terminals.  This  facility  has  been  installed  in  a  separate 
shielded  room  in  our  high  bay  area,  and  is  currently  being  connected  to  the 
various  experiments  it  is  designed  to  support.  In  addition  to  this 
facility,  a  40  kV,  40  A,  5  ys  pulse  modulator  has  been  installed  to  power 
the  CUSPTRON  experiment,  and  this  new  modulator  is  now  fully  operational. 
Additional  remarks  about  the  experimental  facilities  can  be  found  in 
Section  III  of  this  progress  report. 

2.  High  Power  Microwave  Experiments 

Experiments  designed  to  use  the  high  voltage  rotating  beam  facility 
have  produced  significant  results  during  the  past  year.  These  experiments, 
in  which  a  rotating  electron  beam  interacts  with  the  modes  of  a  magnetron- 
type  conducting  boundary  system,  have  demonstrated  that  efficient  radiation 
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can  be  produced  using  either  annular  sector  resonator  or  hole  and  slot 
resonator  configurations.  Using  a  2  MeV,  1-2  kA,  5  ns  rotating  electron 
beam  pulse,  about  300  MW  have  been  produced  at  8  GHz  using  a  hole  and  slot 
resonator  system  with  ten  resonators,  and  about  500  MW  have  been  produced 
at  15.5  GHz  using  a  20  slot  annular  sector  resonator  system.  These  results 
are  detailed  in  a  paper  entitled,  "High  Power  Microwave  Generation  from  a 
Large  Orbit  Gyrotron  with  Vane  and  Hole  and  Slot  Resonators,"  published  in 
the  IEEE  Trans.  Plasma  Science  (special  issue  on  high  power  microwave 
generation)  and  enclosed  in  Appendix  B. 

3.  Millimeter  Wave  Generation  from  Rotating  Beams  in  Rippled 

Magnetic  Fields 

These  experiments,  conducted  in  collaboration  with  Professor  George 
Bekefi's  group  at  MIT,  are  designed  to  investigate  a  novel  kind  of  free 
electron  laser  configuration.  In  this  work,  a  rotating  beam  interacts  with 
an  azimuthally  periodic  wiggler  field  provided  by  samarium  cobalt  magnets 
located  interior  and  exterior  to  the  beam.  Narrow  band  radiation  at  power 
levels  in  excess  of  200  kW  has  been  observed  at  90  GHz  and  180  GHz  in  these 
experiments,  and  the  studies  are  detailed  in  a  paper  entitled,  "Millimeter 
Wave  Generation  from  a  Rotating  Electron  Beam  in  a  Rippled  Magnetic  Field," 
published  in  the  Physics  of  Fluids  and  enclosed  in  Appendix  B. 

4.  CUSPTRON  Experiments 

These  experiments  on  our  low  voltage  table  top  rotating  beam  facility 


(25-40  kV,  1-2  A,  5  ps,  100  pps)  are  now  fully  underway.  Over  500  W  of 
microwave  power  at  the  sixth  harmonic  of  the  cyclotron  frequency  have  been 


observed  when  the  rotating  beam  interacts  with  a  six  slot  vane  resonator 
system,  and  experiments  designed  to  produce  radiation  at  the  twelfth 
harmonic  are  currently  in  progress.  The  potential  of  this  new  device  lies 
in  the  extremely  low  magnetic  fields  required  for  its  operation  (typically 
200  gauss  compared  to  many  kilogauss  for  comparable  magnetrons  or 
gyrotrons) . 

B.  Theoretical  Research 

The  theoretical  contributions  to  our  research  this  past  year  have  been 
the  best  since  the  inception  of  our  program.  Not  only  have  details  of  the 
many  calculations  that  have  been  done  by  our  group  merged  into  common 
agreement,  but  our  ideas  as  to  how  systems  should  be  designed  to  enhance 
radiation  have  been  quantified.  We  believe  that  we  can  now  delineate  the 
main  factors  that  determine  the  best  conditions  for  optimum  coupling  of  a 
rotating  beam  to  a  given  wall  geometry.  The  theoretical  areas  studied 
were : 

•  The  stability  of  a  thin  E  layer  via  a  Vlasov  formalism  was  completed 
and  shown  to  agree  with  the  single  particle  approach.  This  was  used 
to  study  the  interaction  of  a  thin  rotating  E  layer  with  a  vane 
resonator  wall  structure  (see  Refs.  23  and  28  of  Appendix  A). 

•  It  was  proved  that  for  a  cusp  injected  system  the  2 x  mode  of  the 
waveguide  structure  is  the  dominant  mode  of  operation.  Thus,  the 
injection  radius  of  the  beam  selects  the  proper  radial  mode  number  of 


the  "2ir"  harmonic,  and  for  strong  coupling  the  beam  must  be  close  to 
the  slotted  wall  structure.  However,  most  of  the  energy  is  coupled 
into  the  fundamental  mode  (see  Ref.  28  of  Appendix  A). 

Both  resonant  and  nonresonant  interaction  is  considered  via  the 
Vlasov  and  fluid  formalisms.  For  the  nonresonant,  long  wavelength 
limit,  the  traditional  negative  mass  instability  is  recovered  (see 
Refs.  23,  29,  and  30  of  Appendix  A). 

The  derivation  of  a  general  growth  rate  expression  for  systems 
involving  large  orbit,  thin,  annular,  relativistic  electron  beams  was 
completed.  This  growth  rate  formula  is  written  in  terms  of 
parameters  which  contain  the  required  information  about  the  waveguide 
wall  structure.  This  method  extends  previous  results  to  systems 
whose  empty  waveguide  modes  are  hybrids.  As  an  example  of  this 
latter  case,  growth  rates  for  coaxial  systems  with  dielectric  liners 
are  calculated  via  this  formalism.  The  use  of  dielectric  liners  may 
be  very  useful  for  low  energy  systems  where  it  is  difficult  to 
simultaneously  have  resonant  interaction  at  a  strong  coupling 
position  of  the  beam  (see  Ref.  31  in  Appendix  A). 

The  production  of  high  power  microwave  radiation  from  a  large  orbit 
gyrotron  in  azimuthally  periodic  boundary  systems  was  studied 
theoretically.  Linear  growth  rates  are  calculated  for  the  2ir  modes 
of  magnetron-like  vane  resonator  (VR)  and  hole-and-slot  resonator 
(HASR)  systems  using  a  general  growth  rate  formalism.  This  formalism 
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was  used  to  design  the  resonator  structures  for  the  experiments 
performed  in  our  laboratory.  The  radiation  characteristics  are  in 
reasonably  good  agreement  with  the  theoretical  predictions  especially 
with  respect  to  radial  mode  competition  problems  (see  Ref.  32  in 
Appendix  A). 

•  The  macroscopic  fluid  equations  were  used  to  study  the  effect  of  beam 
thickness  on  the  linear  growth  rates  of  large  orbit  gyrotron 
configurations.  Two  methods  were  used  to  calculate  the  growth 
rates.  The  analytic  approach  solves  the  problem  to  first  order  in 
the  ratio  of  beam  thickness  to  the  average  beam  radius  and  recovers 
the  zero  thickness  limit.  The  numeric  approach  uses  complete 
orthonormal  systems  comprised  of  the  homogeneous  wave  equation 
solutions  to  calculate  the  growth  rates  for  beams  with  arbitrary 
thicknesses.  It  is  concluded  that  thick,  tenuous  beams  in  resonant 
systems  must  be  operated  near  the  cutoff  of  the  EM  wave.  The 
decrease  in  growth  rate  as  the  wavelength  decreases  is  attributed  to 
the  loss  of  synchronism  in  systems  where  the  beam  has  a  spread  in  its 
canonical  angular  momentum  [see  Ref.  33  in  Appendix  A  (paper  in  the 
process  of  being  written)]. 

•  The  study  of  single  particle  motion  in  hollow-cylindrical  waveguide 
fields  has  been  initiated.  A  purely  linearized  analytic  approach  can 
be  used  for  particles  far  from  resonance  and  an  analytic-iterative 
analysis  can  be  used  for  particles  that  are  near  or  at  resonant  if 
they  are  trapped.  These  analytic  results  are  compared  to  a  numerical 


solution  of  the  equations  of  motion.  The  main  information  in  which 

we  are  interested  pertains  to  the  process  of  physical  bunching  and  to 

the  net  energy  gained  or  lost  by  a  representative  group  of  particles 

in  the  electromagnetic  field.  We  define  the  magnitude  of  physical 

bunching  by  the  percentage  of  particles  trapped  in  their  azimuthal 

wave  period.  The  effect  of  the  magnitude  and  slope  of  the  electric 

field  is  also  investigated.  For  our  set  of  parameters,  we  determine 

that  dy/dt  is  solely  responsible  for  azimuthal  bunching.  This  can  be 

seen  in  both  the  linear  analytic  theory  and  the  numerical 

solutions.  We  also  observe  that  the  tangential  magnetic  force  causes 

a  spreading  of  the  bunch  in  the  z  direction.  This  spreading  does  not 

occur  for  k  =0  and  v  =0,  i.e.  cutoff,  since  in  that  case  the 
z  zo 

tangential  magnetic  force  is  zero.  Energy  loss  (gain)  by  the 
particles  is  seen  to  be  related  to  both  the  percentage  trapped  and 
the  position  of  the  bunch  in  the  phase  of  the  wave  [see  enclosed 
abstract  (Ref.  35  in  Appendix  A)]. 
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Microwave  radiation  can  be  produced  at  cyclotron  harmonics  as  a  result  of  the  synchronous 
interaction  between  a  rotating  E  layer  and  an  em  wave.  The  derivation  of  a  general  growth  rate 
expression  is  presented  here  for  systems  involving  large  orbit,  thin,  annular,  relativistic  electron 
beams.  This  growth  rate  formula  is  w  ritten  in  terms  of  parameters  which  contain  the  required 
information  about  the  waveguide  wall  structure.  This  method  extends  previous  re  alts  to  systems 
that  have  empty  waveguide  modes  which  are  hybrids.  As  an  example  of  this  latter  case,  growth 
rates  for  coaxial  systems  with  dielectric  liners  are  calculated  via  this  formalism. 


I.  INTRODUCTION 

Microwave  generation  at  cyclotron  harmonics  via  large 
orbit,  annular  electron  beams  has  received  considerable  at¬ 
tention  during  the  past  few  years.  Numerous  theoretical  and 
experimental  efforts  have  resulted  in  a  good  understanding 
of  the  radiation  process.1-12  Studies  have  been  done  for 
beams  in  various  waveguide  geometries  with  the  magne¬ 
tron's  vane  resonator  system  receiving  most  of  the  attention. 

One  promising  device  at  the  University  of  Maryland 
generates  a  rotating  beam  with  a  magnetic  cusp  field.11  This 
device  has  been  used  to  generate  power  levels  of  500  MW 
both  at  the  12th  and  20th  harmonics  of  the  relativistic  cyclo¬ 
tron  frequency  with  a  vane  resonator  system. 1415  The  hollow 
E  layer  nominally  has  a  particle  energy  of  2  MeV,  a  current 
of  1-2  kA,  and  a  pulse  length  of  5-10  nsec.  Other  experi¬ 
ments  on  cusp-injected  devices  are  in  progress  at  Maryland1'1 
and  elsewhere. 1 7  Furthermore,  other  injection  schemes  have 
successfully  been  developed.11* 

Several  of  the  previous  linear  growth  rate  analyses  used 
the  fields  explicitly  associated  with  the  particular  waveguide 
in  the  analysis  (cf.  Refs.  1  and  2).  This  makes  it  necessary  to 
repeat  most  of  the  procedure  every  time  the  waveguide  ge¬ 
ometry  is  changed.  In  fact,  a  large  part  of  the  analysis  is 
independent  of  the  waveguide  and  should  only  have  to  be 
done  once.  In  this  analysis  we  perform  all  of  the  geometry- 
independent  work  first.  Some  geometry  parameters  are  then 
defined  that  enable  us  to  write  a  growth  rate  formula  without 
a  priori  knowledge  of  the  waveguide  wall  structure.  For  a 
given  system  the  geometry  parameters  can  be  found  by  fol¬ 
lowing  a  simple  procedure  involving  only  the  application  of 
boundary  conditions  that  do  not  involve  the  beam  (similar  to 
Ref.  1 1).  Consequently,  we  have  reduced  this  linear  growth 
rate  problem  to  an  algebraic  exercise.  As  an  added  feature 
this  procedure  will  generate  the  dispersion  relation  for  the 
system  without  the  beam. 

Up  to  this  time  all  the  linear  growth  rate  studies  have 
been  valid  only  for  systems  that  have  empty  waveguide 
modes  which  are  TE  or  TM.  In  addition  to  recovering  pre¬ 
vious  results,  our  method  works  on  systems  that  have  empty 
waveguide  modes  which  are  hybrid.  This  enables  us.  for  ex¬ 
ample,  to  calculate  growth  rates  for  waveguides  with  dielec¬ 
tric  liners. 


The  growth  rate  formula  is  derived  in  Sec.  II.  The  geom¬ 
etry  parameters  are  defined  and  the  special  cases  of  trans¬ 
verse  electric  and  magnetic  waves  are  discussed  as  well.  In 
Sec.  Ill  we  recover  the  results  of  Sprangle  (Ref.  1)  and 
Destlerefa/.  (Ref.  2)  and  discuss  the  possible  need  for  dielec¬ 
tric  liners.  Growth  rates  for  a  coaxial  system  with  a  dielectric 
liner  on  the  outer  wall  are  derived  in  Sec.  IV.  The  special 
cases  for  TE  and  TM  waves  are  carefully  discussed  before 
the  general  case  is  derived.  Conclusions  are  drawn  in  Sec.  V. 

II.  THE  GENERAL  GROWTH  RATE  FORMULA 

The  linear  growth  rate  problem  can  roughly  be  divided 
into  two  parts.  The  first  part  involves  calculating  the  per¬ 
turbed  sources  in  terms  of  the  em  fields  at  the  equilibrium 
radius.  Because  this  calculation  has  essentially  been  done 
before  (cf.  Ref.  1),  we  only  briefly  sketch  the  derivation  here. 
In  this  sketch  we  reorganize  the  equations  to  better  fit  the 
general  analysis  and  emphasize  the  key  points  in  the  deriva¬ 
tion. 

In  the  second  part  we  close  the  system  of  equations  by 
integrating  two  components  of  Ampere’s  law  across  the 
beam.  This  produces  two  equations  that  express  thejumps  in 
cB,  and  r  d,  Et  across  the  beam  in  terms  of  certain  beam 
parameters  and  the  values  of  rdrcB ,  and  E,  at  the  beam 
radius.  It  is  at  this  point  that  the  waveguide  geometry  first 
enters  into  the  analysis.  The  particular  boundary  conditions 
at  the  waveguide  walls  determine  the  ratios  of  cB,  and 
r drE,  to  r d.cB,  and  E,  at  the  beam  radius.  In  fact,  the 
geometry  parameters  used  in  the  growth  rate  expression  are 
composed  of  these  ratios.  These  geometry  parameters  are 
generalizations  of  the  (two)  normalized  wave  admittances 
used  by  others  in  the  analysis  of  decoupled  systems  (see  Ref. 
3,  for  example).  The  generalization  is  necessary  because  the 
usual  wave  admittances  are  insufficient  for  the  hybrid  case. 
The  reformulation  of  the  em  fields  in  the  source-free  regions 
and  the  incorporation  of  the  geometry  parameters  are  the 
main  new  features  in  (his  analysis. 

A.  The  perturbed  source  terms 

The  type  of  system  we  consider  has  an  electron  beam 
propagating  in  some  cylindrical  waveguide  which  is  im- 


2868  Ptiys.  Fluids  28  (9),  September  1985 


0031-9171/85/092868-10S01  90 


©  1985  American  institute  ot  Pttytcs 


2868 


FIG.  1.  The  equilibrium  model. 


mersed  in  a  uniform  magnetic  field  (see  Fig.  1).  The  only 
restrictions  placed  on  the  waveguide  are  that  it  must  be  axi- 
symmetnc  and  invariant  to  translations  in  the  z  direction. 
These  restrictions,  along  with  the  linearization,  allow  us  to 
assume  that  only  one  harmonic  perturbation  exists  in  the 
system.  That  is,  perturbations  with  different  harmonic 
numbers  do  not  interact  in  any  way. 

We  assume  that  the  beam  is  very  thin;  its  equilibrium 
density  is  given  by 

n0{r)  =  n,S(r  -  r0), 

where  r0  is  the  beam  radius  and  n,  is  the  surface  particle 
density.  We  also  assume  that  the  beam  is  sufficiently  tenuous 
so  that  the  equilibrium  self-fields  can  be  neglected.  Each 
electron  moves  on  a  predominantly  helical  orbit  with  an  un¬ 
perturbed,  normalized  velocity  of  (0,  Pe  ,  /?4  ).  This  simple 
distribution  enables  us  to  use  the  single  particle  equations  in 
our  analysis. 

The  first  step  is  to  linearize  and  Fourier  analyze  the  Lor- 
entz  force  equation.  We  write  the  particle  positions  as 

r  =  r0  +  r„  9  =  90  +  fl0t  +  6„  z  =  z0  +  vj  +  z„ 

where  /?„  =  eB,/{mn  y0)  is  the  electron-cyclotron  frequency 
and  Y(/n0 cJ  is  the  electron  energy.  The  perturbations  r„  9„ 
and  z,  are  assumed  to  be  small  compared  to  their  zeroth- 
order  counterparts. 

We  assume  that  all  perturbations  have  their  dependence 
on  the  time,  f,  and  the  9  and  z  coordinates  defined  by 
exp[/(fc,z  + 19  —  at)].  The  amplitudes  of  the  perturbations 
are  denoted  6„  and  z,.  Consequently,  the  linearized, 
Fourier  analyzed  version  of  the  total  time  derivative  is 
d  /dt—  —  til/, ,  where 

tft,=a-in0-k,v^.  (1) 

The  quantity  i>,  is  a  measure  of  the  difference  between  the 
actual  em  wave  in  the  beam-waveguide  system  and  the 
doppler-shifted  cyclotron  harmonic.  In  the  following  analy¬ 
sis  lit,  plays  a  very  important  role  .  The  synchronous  limit  is 
defined  by  tf,— *0.  Thus  we  say  that  the  beam  is  synchronous 
with  the  em  wave  when  \l>t  is  small  compared  to  the  other 
frequencies  in  the  analysis  (if',  </?„  in  particular).  In  our  deri¬ 
vation  we  always  assume  the  synchronous  limit  and  only 
keep  the  lowest  order  terms  in  t .  This  assumption  has  many 
consequences.  Most  importantly,  it  enables  us  to  obtain  an 
explicit  growth  rate  formula.  The  other  consequences  will  be 
discussed  as  they  appear.  We  will  show  later  that  the  syn¬ 
chronous  assumption  is  a  good  one  for  tenuous  beams. 


The  linearized,  Fourier-analyzed  Lorentz  force  equa¬ 
tion  is  found  to  be 

r,  =  (  -  ie/m0  yji^,  )(£„  +  P,ocB, )  +  O  ( 1 ),  (2a) 

ro0,  =  (  -  e/m0  y0\/>f)(p2oEe  +  P6u  p^E,)  +  0(^y  '),  (2b) 
z,  =  (e/m0  y0ttf)  [  ( 1  ~P]JE,  -  PecB,  -  P6o  P^E#].  (2c) 

Equation  (2)  reveals  that,  to  lowest  order  in  ilr, ,  the  perturbed 
particle  positions  depend  only  on  the  field  components, 
which  are  continuous  across  the  beam.  This  is  convenient 
because  the  continuous  fields  are  well  defined  at  the  beam 
radius  in  terms  of  limiting  values  as  r— *r0.  When  r^r0,  the 
vacuum  equations  apply  and  the  continuous  fields  can  be 
written  in  terms  of  E,  and  r  3,cB, . 

Conservation  of  electrons  on  each  differential  surface 
element  requires  that,  again  to  lowest  order  in  i/>,  (cf.  Ref.  1 ), 

n,  =  -  in, (10,  +  k,z,)S(r  -  r0).  (3) 

From  this  we  find  the  perturbed  current  to  be  (as  always,  to 
lowest  order  in  tl>, ) 

#y,,/f(lc  =  0,  (4a) 

oC=  -  {r(&b/c)2f}(l 

X  [  -  E  ;  +  n0(r  3 ,  cB , )°  ]  6(r  -  r0). 

(4b) 


X  [  -  iI2mE°,  +  nu(rd,cB,)°]8(r  -  /•„), 

(4c) 

where  a\  -  e2ns/(m0y0r„e0),  /?,  =(ktc  —  ap^),  and  the 
superscript  zeros  mean  that  the  function  is  evaluated  at  the 
equilibrium  radius  r0.  From  Eqs.  (3)  and  (4)  we  see  that,  in  the 
synchronous  limit,  the  source  perturbations  lie  on  the  r  =  r0 
surface.  This  fact  will  be  important  in  the  next  section  when 
we  apply  the  jump  conditions. 


B.  The  jump  conditions  and  the  geometry  parameters 

There  are  two  unknowns  on  the  rhs  of  Eq.  (4)  and  there¬ 
fore  two  more  equations  are  needed  to  close  the  system.  We 
integrate  the  Fourier-analyzed  9  and  z  components  of  Amp¬ 
ere’s  law  across  the  beam  to  get 

-  75  _  P*  Jo.  \r)dr 


5.£- r"-- —■ 

0  Jr„-  rn  fnC 


"  J'o  f<)  foC 

Substituting  Eq.  (4)  into  Eq.  (5)  and  rewriting  cBe  in  terms  of 
the  axial  fields  gives 

tfMr£,)|J  =2-^(  — ,/2.) 

Vo 

X  [  —  if2mE°  +  fl0(r 3,cBt)n)y  (6a) 

*?(<*,)£  =  2  - («„)[  - 1/2, E°,  +  ajrd'CB, )°J. 

To 


where  v  =  r0e2n,  /i,/(2m0)  is  Budkcr's  parameter.  Nothing 
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more  can  be  done  without  considering  the  waveguide  geome¬ 
tries. 

When  r#r0  the  vacuum  waveguide  equations  apply.  In 
cylindrical  coordinates,  the  general  solution  is  a  linear  com¬ 
bination  of  Bessel  functions  of  the  first  ( J , )  and  second  ( >', ) 
kinds.  We  divide  the  r  direction  into  two  regions:  region  I, 
where  r<r0,  and  region  II,  where  r>r0.  Symbolically,  we 
write  the  axial  fields  in  the  two  regions  as 


J.  _  |/„  ('•)£“  +  '/,»  ('■)('■  d,cBr )°,  r  <  r„ 
\fZ[r)E°,  +  ‘f.l  (''H'’  drCBzf,  r>ru, 


(7) 


and 


'/ b.  ('■)£“  +/»  ('■)('•  d,cB, )°,  r  <  r0, 
if  A  ('')£“  +/m  ('•)('■  drCB,  )u.  r  >  r0. 


The  radial  functions //,  (r)  are  just  linear  combinations  of  Jt 
and  Y,.  Their  coefficients  are  determined  by  boundary  con¬ 
ditions  that  do  not  involve  the  beam.  Examples  of  the  calcu¬ 
lations  of  these  functions  are  given  in  the  following  two  sec¬ 
tions. 

To  proceed  with  the  general  calculation,  we  define 


<*4 

—  f  bb  l^o)  ~fbb  ( 

t,  =  0 

m  | 

ird,cB,r 

ro) 

_  '**4 

dn 

df~ 

C  •, 

1 

IrO.cB.f  =  O  dr 

ro 

dr 

<*4 


St* 


‘E°, 


—  f  bt  (ro)  ~fbt  (ro)> 


=  0 


Seb  .. 


rd,EA% 

dffb 

—  jr 

df,b 

i(rd,cBzf 

£?-  o  dr 

'«  dr 

(9a) 


(9b) 


(9c) 


(9d) 


to  be  the  geometry  parameters.  We  label  gbb  the  TE  param¬ 
eter,  g„  the  TM  parameter,  and  and  gtb  the  hybrid  pa¬ 
rameters.  Plugging  Eq.  (9)  into  Eq.  (6)  yields  the  dispersion 
relation 


#(&,!,  g„  +  gb,  grh  ) 

=  2{v/r0)[  g„ni  -  (g*  +gtb)n,n0 -gbbn\ ] 

(10) 

for  the  beam  waveguide  system  in  the  synchronous,  tenuous 
beam  limit. 

When  there  is  no  beam,  v  =  0,  and  the  empty  waveguide 
dispersion  is 

D.(oj,k, )  =  gbb  g„  +  gh,  g,b  =  0.  (11) 

We  denote  solutions  of  Eq.  (1 1)  by  (to*  ,k't).  If,  by  a  suitable 
choice  of  system  parameters,  the  perturbed  beam  frequency 
lfin  -f  k’v^  is  equal  to  some  empty  waveguide  frequency 
to’ ,  we  say  that  the  beam  perturbation  is  in  resonance  with 
the  empty  waveguide  mode  (at  k ').  Resonant  interaction  for 
a  simple  waveguide  geometry  is  depicted  in  Fig.  2.  Wc  note 
that,  in  general,  there  are  two  resonant  points,  if  any.  The 
intersection  corresponding  to  k'  t[k’_  )  is  called  the  high- 
flow-)  frequency  intersection  point.  If  k  J  _  <0,  then  the  res¬ 
onance  is  called  a  backward  wave  intersection;  otherw  ise  it  is 


FIG.  2.  Graphical  description  of  resonant  interaction. 


a  forward  wave  intersection  [k‘  =  0  is  called  the  cutoff).  For 
a  tangential  intersection,  there  is  only  one  resonant  point  and 
that  point  is  characterized  by  identical  group  velocities  for 
the  waveguide  and  beam  modes,  i.e.,  the  grazing  condition. 

C.  Properties  ot  the  growth  rate  formula 

To  keep  the  analyi  s  consistent,  the  geometry  param¬ 
eters  must  also  be  evaluated  to  lowest  order  in  ip, .  Although 
Dt((o,kz )  is  usually  zeroth  order  in  1 1{,  at  resonance  it  be¬ 
comes  first  order.  Therefore,  at  resonance  the  dispersion  re¬ 
lation  becomes  a  cubic  in  instead  of  a  quadratic.  In  this 
paper  we  consider  only  the  resonant  case  because  the  cubic 
growth  rate  dominates  the  quadratic  growth  rate  at  low  den¬ 
sities. 

We  emphasize  that  the  dispersion  function  [Eq.  (10)]  is 
valid  off  resonance.  Only  synchronism  is  required  by  the 
analysis.  For  example,  in  the  appropriate  limit  we  could 
show  that  our  growth  rate  agrees  with  Uhm's  negative  mass 
growth  rate  in  coaxial  geometry.19 

At  resonance  we  obtain  from  Eqs.  (10)  and  (11)  that 

=  2  fil)  g”n  °  ~ ( 8bf  +  S'tf2"0*  ~  • 

W  (£<*£„  +gb'gtbY 

(12) 

where  the  prime  denotes  the  derivative  with  respect  to  to  and 
the  geometry  parameters  have  been  evaluated  at  [to’  ,k ').  We 
note  that  the  frequency  shift  tf>,  depends  only  on  v/y„,  the 
frequencies  and  ,  and  the  geometry  parameters.  We 
write  il’i  —  to,  +  iT  so  that  /">  0  represents  a  growth  in 
beam  amplitude  with  time.  Assuming  the  rhs  of  Eq.  (12)  is 
real,  the  growth  rate  and  frequency  shift  are  given  by 

r  =  (4/ 2)|<M,,\  (13a) 

to,  =  —  sgn(tf,')/~  A/3,  (13b) 

where  sgn  denotes  the  signum  function. 

It  follows  from  Eqs.  (6)  and  (9)  that  the  empty  waveguide 
modes  are  TE  or  TM  if gv  =  0  =  grb .  In  this  case  the  empty 
dispersion  relation  and  corresponding  growth  rate  for  a  I  E 
mode  arc 
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**=  0.  (14a) 

£**  =  ,14b) 

o*  2  \a0eJ 

The  corresponding  expressions  for  a  TM  mode  are 

g_  =0.  (15a) 

r^r.  (15b) 

a.  2  \a.fj 


Notice  that  the  TM  growth  rate  depends  on  J?c,  but  the  TE 
growth  rate  depends  on  flo-  Consequently,  iffl,  =0,  then 
the  TM  modes  are  stable  but  the  TE  modes  are  still  unstable. 
This  can  be  seen  physically  by  considering  a  transformation 
to  the  beam  frame  (moving  with  a  velocity  with  respect  to 
the  lab  frame).  For  a  TM  mode/?,  =  0  implies  that  EL  =0 
in  the  beam  frame:  Thus  y  is  second  order  and  particles  do 
not  lose  energy  in  the  linear  approximation.  For  a  TE  mode, 
however,  Et  is  finite  in  the  beam  frame  and  y  remains  a  first - 
order  term.  The  simplest  way  to  obtain  12 0  =  0  is  to  set 
a^  =0  =  4,.  In  the  dielectric-free  cases,  fl,  =0  whenever 
the  grazing  condition  is  satisfied 

Also,  notice  that  /'<x  n'/1.  This  scaling  law  implies  that 
the  synchronous  limit  is  valid  for  sufficiently  tenuous  beams. 
In  the  next  section  we  use  Eqs.  (14)  and  (15)  to  quickly  re¬ 
cover  the  results  of  Sprangle  and  Destler  el  al. 


FIG.  3.  Cross  section  of  the  simple  cylindrical  waveguide  geometry. 


and  j 

/«i  =/rt  =/£  . 

%> 

for  f  2  =  {(it/c)2  —  k),  x(r)  =  r£,  and  x0  =  x(r0),  etc.  From 
Eq.  (9)  the  TE  and  TM  geometry  parameters  are  found  to  be 

_  2  •*?(*,) 
gU  *  xS/;(x0)S/(*o**J 

and 


OL  DISCUSSION  OF  PREVIOUS  WORK 

In  this  and  the  next  section  it  is  convenient  to  adopt  the 


standard  notation  for  the  cross  products  of  Bessel  func- 

tinrtv-20 

7 

Pfcj> )  =  JMYi(b )  -  Y,{a\I,(b ), 

(16a) 

QM* )  =  J,l*)r;\b )  -  Y, (a)/;(ft ). 

(16b) 

R,{pJb)^J](a)Y,(b)- T,W,(H 

(16c) 

S,1pjb )  =  J;(a)Y](b )  -  Y'AaU \{b ). 

(16d) 

_  _  2  -M*J 

g"  w  JAx0]P,{Xo^j' 

As  expected,  the  TE  solutions  are  characterized  by  J  «'(xw) 
=  0  and  the  TM  solutions  are  characterized  by  J,[x„)  —  0. 
The  growth  rates  are  found  from  Eqs.  (14)  and  (15)  to  be 

r te  =  ^  L  v  (c/rj^y^ 
no  2  V  y0  (l&>  ) 

I _ *«/«(*o) _ I171 

I[i-(/a.)2],/v,(xm)| 


It  is  simple  to  show  from  Bessel’s  equation  that,  if  a  =  a{z) 
and  b  =  b  (z),  then 

P]\flJb)  =  aR,  +b'Q„  (17a) 

C;M)  =  fl'5f  +*'((/ Vh 2 -  1JP,  ~(\/b)Q,\.  (17b) 
Jt;(fl^)  =  a*[(/2/o2-  1JF,  ~(\/a]R,\  +b'S„  (17c) 

)  =  «'[(/ Va2  -  1  )&  -  (l/a)S,  ] 

+  h'[(/2/62-l)A,-(l/h)S/].  (17d) 

Clever  choices  for  arguments  of  these  cross  products  sim¬ 
plify  the  algebra  tremendously. 

Consider  the  case  of  only  a  smooth  outer  conducting 
wall  of  radius  rm  (see  Fig.  3).  We  apply  the  usual  boundary 
conditions  at  rm  and  growth  conditions  at  r  =  0  and  make 
the  following  identifications  (via  Eqs.  (7)  and  (8)]: 

JM  ,  >  pfc<x.) 

im  J.w’  Jm  P.^xJ  ’ 

_  JM 

/m  “  ~  r 9  j  *  ~  zn  r  *  * 

VIM 


and 

£Vm  =  I  J'M  | 2/2 

2  \  y0  I2ma>  )  |  J](xJ  | 

These  results  do  indeed  agree  with  Sprangle’s  synchronous 
case  (Ref.  1). 

When  a  smooth  inner  conductor  of  radius  r,  is  added  to 
the  system  (see  Fig.  4),  only  the  boundary  conditions  in  the 
region  below  the  beam  are  affected.  The  net  difference  is  in 
the  two  functions  f  „  and  f  ^ ,  which  become 

fm  ~  P t^XiXiVP f(*o**/) 

and 

fu  =  Ql  {X*i  )/XaS,  (x^, ) 

for  the  inner  conductor  case.  The  geometry  parameters  be¬ 
come 

g  -  2  ^(x.^,) 

*  v  x$S(  (*„,*,  WwJ 

and 
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FIG.  4.  Cross  section  of  the  coaxial  waveguide  geometry. 


_  _  2  Pt(xwjc,) 

8“  *  /'ilvJ/’ilvJ 

so  that  TE  solutions  are  characterized  by  S,{xu )  =  0  and 
TM  solutions  by  P,{xwjc,)  =  0.  The  TE  and  TM  growth 
rates  are  found  from  Eqs.  (14)  and  (IS)  to  be 


A, 

£tm 


V3  l_v  (^j2y/} 

2  V  y0  ) 


{rr/Dx^.x'oS,  [xajc,  )5,  (x„jcw ) 
dS,{xwrx,)/dxw 


t/j 


X 


{ir/2pcwP,{x0Tx,  )P,{x„,xu) 
dP,(xwj t,  )/dxw 


i/j 


where  Eq.  (17)  is  used  to  compute  the  derivatives. 

These  results  agree  with  Destler  (Ref.  2)  except  for  an 
amplitude  factor.  As  pointed  out  by  Destler  et  ai,  the  dis¬ 
crepancy  comes  from  the  method  of  calculating  the  per¬ 
turbed  charge  density.  There  is  strong  evidence  pointing  to 
the  correctness  of  our  results.  First,  it  was  already  pointed 
out  here  that  our  results  agree  with  Uhm's  kinetic  theory 
approach  in  the  proper  limit.  Furthermore,  a  fluid  analysis 
done  by  our  group:i  has  recovered  the  dispersion  relation 
Eq.  (10)  in  the  infinitesimal  E  layer  limit  (at  least  when  g 

=  0  =  g,b ). 

In  both  the  cases  with  and  without  an  inner  conductor, 
if  resonance  is  possible  then  the  resonant  wavenumbers  are 
given  by 

I ck,  )t  =ri ma P„  ±  Y,„  [  (n„ M0)'  -  (i-c)1] 1/2  (18) 

for  1  =  1  —  p  j .  The  condition  for  resonance  can  be  writ¬ 
ten 


r„  l  YoPn, 

rw  <  xw  (1  +  (YoPo.,)1]'1 


(19) 


At  high  azimuthal  harmonics,  the  peak  Ee  field  is  very  close 
to  the  wall.  Thus  in  order  to  get  reasonable  growth  rates,  it  is 
necessary  to  make  the  ratio  r,/r,u  as  large  as  possible.  For  all 
possible  TE  or  TM  solutions  of  either  case,  it  is  always  true 
that  I  <xw.  Since  P„a  <  \,r0/rw  must  be  very  small  unless  the 


beam  energy  is  large.  Consequently,  we  do  not  expect  low- 
energy  systems  to  be  very  successful  in  these  simple  geome¬ 
tries. 

For  low-energy  systems  something  must  be  done  to  (in 
essence)  raise  the  possible  values  of  /  /xw .  This  is  just  another 
way  of  saying  that  something  is  needed  to  “slow  down"  the 
em  wave.  One  possibility  is  the  slotted  wall  structure.  This 
method  has  received  considerable  attention  and  will  not  be 
discussed  here.  Another  possible  method  involves  the  inser¬ 
tion  of  a  dielectric  liner.  This  approach  is  studied  in  detail  in 
the  next  section. 


IV.  THE  DIELECTRIC  LINER 

Consider  the  system  shown  in  Fig.  5.  It  is  similar  to  the 
coaxial  system  except  that  a  dielectric  liner  of  relative  per¬ 
mittivity  e ,  occupies  the  space  between  r  =  r4  and  r  =  rw. 
As  before,  the  beam  radius  is  r0  and  the  inner  and  outer  wall 
radii  are  r,  and  rw,  respectively.  There  is  no  fundamental 
difference  between  this  system  and  one  without  an  inner  con¬ 
ductor.  In  fact,  the  results  are  virtually  identical  until  the 
inner  radius  becomes  relatively  large.  However,  the  disper¬ 
sion  relation  for  the  coaxial  dielectric  liner  case  takes  on  a 
more  symmetric  form  and  is  easier  to  solve  numerically.  For 
these  reasons  we  present  the  details  of  the  calculation  for  the 
coaxial  case.  Afterward,  we  present  the  results  for  the  special 
case  r,  =  0. 

The  procedure  for  obtaining  the  geometry  parameters  is 
still  straightforward,  but  now  the  algebra  is  more  complicat¬ 
ed.  Also,  the  results  are  harder  to  analyze.  First,  the  modes 
only  decouple  into  TE  and  TM  modes  at  isolated  points. 
Second,  the  dispersion  curves  are  no  longer  simple  hyperbo¬ 
las  in  co  —  k.  space.  Finally,  the  arguments  of  some  of  the 
Bessel  functions  become  imaginary  in  part  of  the  a  —  kt 
plane. 

In  Sec.  IV  A,  we  derive  the  geometry  parameters  and 
the  empty  waveguide  dispersion  relation.  In  Sec.  IV  B,  we 
discuss  in  detail  the  special  cases  where  the  modes  are  either 
TE  or  TM.  We  present  the  general  result  and  apply  it  to  both 
a  high-energy  and  a  low-energy  system  in  Sec.  IV  C. 


FIG  5.  Crow  section  of  a  coaxial  system  with  a  dielectric  liner. 
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A.  The  empty  waveguide  dispersion  relation 

In  anticipation  of  the  results  to  come,  we  define 

=(a)/c)2-k],  (20a) 

H  =f,(a>/c)2  -k\,  (20b) 

and 

*'=nok,/c){er-  IV£1.  (21) 

We  also  define  five  sets  of  cross  products: 

=  <£/(■*.*, ).  4>\  =  ^M),  (22a) 

=  ft  Md ).  ^2  =  <^('0).  (22b) 

4>M  =  </>,  (x,x0),  ^3  =  &M ).  (22c> 

$4  =  <*>,(Wi),  (22d) 

=  <*>/  ( y,yw ).  ^5  =  <*,s('d ).  (22e) 

for  <P  =  P,  Q,  R,  and  S;  for  x  —  r£K,  y  =  r£B,  and 
=  rdgB,  etc.  Finally,  we  define  some  more  general  cross 
products: 

Ki-J  =  xd{er;PlR1  -PjR,),  (23a) 

XL'  =  *MQ.%-QjSX  (23b) 

=  V'PQ,,  (23c) 

=  {23d) 

f°r  £  =  ^*/>b- 

Now  we  must  divide  the  r  coordinate  space  into  three 
regions.  We  write  the  axial  fields  symbolically  in  these  re¬ 
gions  as 


’frAr)E°'+ift-bnrdrcB,)°, 

r,<r<r0. 

f:\r)E°  +  i/:b(r)(rd,cB')0, 

r0<r<rd, 

(24) 

f:.\r)E°  +  if%[r)[rdrcB,f, 

rd  <r<,rw. 

VwM®?  +fbb{r){rd,cBzf, 

r,<r<r0, 

•  if£(r\E°  +f:b(r)(rd,cBI)°, 

r0  <r<rd. 

(25) 

ifl{r)E0t  +  /  *j,  (r)(r  <?,  cBz )°, 

rd  <r<rw. 

We  then  systematically  apply  all  the  boundary  condi¬ 
tions  that  do  not  involve  the  beam  in  order  to  identify  the  12 
functions  defined  above.  In  addition  to  the  usual  boundary 
conditions  at  r, ,  r0,  and  rw ,  there  are  six  boundary  conditions 
at  rd .  The  interesting  boundary  condition  is 

E,\t_  =  erE,\r.\ 

the  other  five  reflect  the  continuous  nature  of  the  remaining 
field  components.  Two  of  the  boundary  conditions  are  re¬ 
dundant.  The  remaining  four  enable  the  task  of  finding  the 
12  functions  to  be  completed.  The  results  are 

/*(')= 0,  /*M=0, 

/«  M  =  PtWPf 

/,:  W  =  -  PM  /£  (r)=  -  GiW- 

/:  w  =  l  (-  -  /») . 

/»;  w  =  -4-  (-  aw  -  /») .  t 


/aw=  /tw=  -  -4rG5(r), 

ir  x0i7  •  jt  r?3-5 

2  yJ-5  2 

/iw  =  pM-  /iM  =  -  ~n  &M- 

7T  7/  IT  Xq7] 

Given  the  field  functions,  the  geometry  parameters  are 
found  from  Eq.  (9)  to  be 

*._(lY-&£l  +  -&._4.  aw 

•~J  4Q,n“  *„<?, 

+  (26b) 

\irJ  pyVy>  ps  /», 

**-«-  -(7)'^  |2W 

We  note  that  the  dielectric  empty  waveguide  modes  decou¬ 
ple  into  TE  and  TM  modes  if  r5  5  =  0. 

Two  equivalent  forms  of  the  empty  waveguide  disper¬ 
sion  relation  will  be  used  in  this  section.  The  first  comes  from 
applying  the  definition  [Eq.  ( 1 1 )] 

Vq  Si  v  tr /  x  0j77  l  Fj  />, 

The  numerator  of  this  relation  appears  to  depend  on  the 
beam  radius  r0.  It  should  not,  because  r0  has  no  physical 
significance  when  there  is  no  beam  and  thus  cannot  affect  the 
zeros  of  the  dispersion  relation.  After  more  algebra,  the  dis¬ 
persion  relation  can  be  written  so  that  the  numerator  is  inde¬ 
pendent  of  r0: 

D,(<o,kt)  =  - ~  =0-  ,  ••  (27) 

(ir/2)2XoPlS,r),s 

Thus  empty  waveguide  solutions  are  characterized  by 

ij*s  =  0. 

A  typical  dispersion  plot  is  shown  in  Fig.  6.  The  cutoff 
line  is  defined  by  k ,  —  0.  The  vacuum  and  dielectric  light 


Varuum  Light  Line.. 
Region  A 


Dielectm  Light  Line 


no  too  200  300  400  son 

c*./(2.)  (Gilt) 

FIG.  6.  Typical  dispersion  plot  for  a  dielcctnc  system.  Parameters.  /  =  7, 
r,  =  0.055  m,  rt  =  0.065  m,  r„  =  0.075  m,  and  t,  =  6.0. 
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FIG.  7.  The  effect  of  liner  thickness  on  the  cutolffrequency  of  the  first  radial 
mode.  Parameters:  /  =  7,  r,  =  0.025  m,  rw  =  0.075  m.  The  measure  of  liner 
thickness  is  y  =  \tw  —  r4  )/|r„  -  r, ). 


lines  are  defined  by  ai/{ckz)  =  1  and  a)/{ck, )  =  e~  '/2,  re¬ 
spectively.  We  denote  the  region  in  co  —  k,  space  between 
the  cutoff  line  and  the  vacuum  light  line  as  region  A.  The 
region  between  the  vacuum  and  dielectric  light  lines  is  de¬ 
noted  region  B.  In  region  A,  £A  and  £B  are  both  real.  In 
region  B,  is  still  real  but  5  A  is  purely  imaginary  and  <P , 
through  <PA  are  given  in  terms  of  modified  Bessel  functions. 
Fortunately,  the  geometry  parameters  are  still  real,  so  the 
growth  rate  equation  [Eq.  (13)]  is  still  valid  for  resonant 
points  in  region  B. 

The  dielectric  has  two  major  effects  on  the  dispersion 
curves.  First,  the  cutoff  points  are  lowered  (in  frequency). 
Figure  7  shows  the  relative  cutoff  frequency  for  the  lowest 
dispersion  curve  versus  the  relative  permittivity  for  several 
liner  thicknesses.  The  figure  reveals  that  the  relative  cutoff 
rapidly  approaches  e~ 1/2  as  the  thickness  increases.  Second, 
the  dispersion  curves  are  asymptotic  to  the  dielectric  light 
line.  This  feature  is  more  important  for  linearly  streaming 
beams  (/  =  0)  where  the  beam  line  is  always  below  the  vacu¬ 
um  light  line.  Depression  of  the  cutoff  frequency  is  more 
important  for  the  cases  of  interest  in  this  paper  (/  ^  0). 


B.  The  special  cases  of  TE  and  TM  modes 

From  the  definition  of  r/J  [Eq.  (23)],  we  see  that 
the  modes  will  be  TE  or  TM  if  lPe  =  0,  Q,  =  0,  or  Ps  =  0. 
Equation  (21)  reveals  that  there  are  three  ways  to  have 
lPe  =  0.  If  e,  =  1,  we  are  obviously  back  to  the  coaxial  case 
without  a  dielectric.  The  second  possibility  is  /  =  0.  The  the¬ 
ory  for  this  case  has  been  considered  elsewhere  and  success¬ 
ful  experiments  have  been  performed  by  J.  Walsh  el  a!.22  at 
Dartmouth. 

The  final  possibility  is  k,  =  0.  At  cutoff,  £H  =  s/f7fA 
and  tj*  5  =  k*ax  The  TE  result  is  specified  by  gth  =  0, 
which  is  equivalent  to  y4J  =0.  As  in  the  case  without  a 
dielectric,  only  cBt,  E,,  and  Ee  are  nonzero  at  this  cutoff 
point.  The  TM  result  is  given  either  by  g„  =0otk4>  =0.  At 


the  TM  cutoff  point,  only  Ez,cBr,  and  cB9  are  nonzero. 
At  a  point  where  Q$  =  0, 


Sbb  ~ 


-  (2/yfo 

xoQj£i 


Thus  the  resulting  TE  mode  is  characterized  by 

(?<(*<,)  =  0  =  {?5(yj.  (28) 

This  TE  mode  has  cBz[rd)  —  0.  Consequently,  E, [rd )  =  0 
and  so  Er  is  continuous  at  rd .  The  growth  rate  is  found  from 
Eq.  (14): 


^TE 

_  V3  / 

>  (c/r0)2 ' 

\,/3  {ir/2]xlQ^S, 

n0 

2  \ 

•  Yo  ‘ 

1  dQt/dx0 

When  P5  =  0, 

g„  =  (2/Trjh/hP, 


and  a  TM  mode  has 

?4(xo)  =  0  =  P5(y„).  (29) 

The  TM  mode  is  characterized  by  E,[rd)  =  0,  almost  as 
though  there  were  a  conducting  wall  at  rd .  The  growth  rate 
is  given  by 


r tm 

a. 


J3  (zv  \c/r0)2y> 

2  V  Yo  ft*™  ) 


(jrW*o P& 
dPJdx9 


1/3 


It  is  known23  that  the  zeros  of  /Mxq)  are  all  real,  so  these 
TM  modes  can  only  exist  in  region  A.  Also,  we  can  show  that 
(?4(;t0)  has  no  purely  imaginary  zeros  so  the  TE  modes  all  lie 
in  region  A  as  well.  A  dispersion  plot  highlighting  the  TE 
and  TM  points  is  shown  in  Fig.  8.  Let  (£A  ,£b)  represent  any 
solution  of  either  Eq.  (28)  or  Eq.  (29).  If  £  B  >  A ,  then  the 

pair  (£a,£b)  corresponds  to  the  solution  of  the  dispersion 
relation  that  has 


o>  =  c[(£2B  -£!)/(*, -1  )]'« 

and 

ck.=c[£l  -*,£!)/(*, -1)],/2. 


FIG  8  A  dispersion  plot  revealing  the  locations  of  the  TE  and  TM  points. 
Parameters:  /=  7,  r,  =  0  020  m.  rt  =0  0o5m,  rm  =  0  075  m,  and  e,  =*  4.0. 
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Harmonic  Number  I 


FIG.  9.  The  linear  growth  rate  curves  for  a  high-energy  beam  in  the  simple 
waveguide  geometry:  (a)  the  lowest  TE  mode.  |b|  the  lowest  TM  mode.  The 
thin  line  represents  the  low-frequency  intersection  and  the  thick  line  repre¬ 
sents  the  high-frequency  intersection.  Parameters:  n,  =  lO14  m  ~ 
Pa,  =  0  935, =  0.305,  r0  =  0.06  m,  and  r„  =  0.075  m. 


Furthermore,  if  £B  represents  any  other  solution  [of 
Fj(  yw )  =  0  or  Q>{  )  =  0,  whichever  is  appropriate]  and 

£  B  >£  b-  then  (£A.£B)  also  corresponds  to  a  solution  of  the 
dispersion  relation. 


C.  The  general  result 

The  calculation  of  the  derivative  of  the  empty  wave¬ 
guide  dispersion  relation  is  straightforward  but  tedious.  The 
final  expression  is  quite  lengthy  and  is  not  presented  here. 
For  our  purposes  it  is  sufficient  to  write 

—  ( 

( g„t  g„  +  g*  g,*)  =  ,  ~  ,  ,  ■  ■  .  (30| 

[ir/2yx;,P ,5,7  -5 

when  74,5  =  0.  Once  again  the  prime  denotes  the  derivative 
with  respect  to  ru. 

The  solution  for  the  case  of  r,  =0  can  be  found  by  fol¬ 
lowing  the  same  procedure  outlined  above.  The  solution  can 


■  V  —  V  -  V  —  V  —  w  -  .  »  •  •  _  ■ 


be  written  in  the  same  form,  provided  that  the  following 
substitutions  are  made: 

P \iQi  *'//(.*o)> 

P*iQ*~*J i(xd\<  ’JAxj !■ 

In  the  following,  we  present  results  both  for  r,  =  0  and 

r,  >0. 

The  first  example  is  for  a  high-energy  beam.  The  param¬ 
eters  for  this  example  are/3w  =  0.935,  y„  =  5.46,  r,  =  0,  and 
rc/rw  —  0.8.  A  plot  of  growth  rate  versus  harmonic  number 
is  shown  in  Fig.  9  for  the  lowest  two  dispersion  curves  when 
there  is  no  liner.Without  a  dielectric  liner  there  is  resonant 
interaction  starting  with  /  =  7. 

Figure  10  shows  the  growth  rates  for  the  system  when  a 
1  cm  thick  dielectric  liner  with  a  relative  permittivity  of  2.3  is 
placed  against  the  outer  wall.  Although  nothing  drastic  hap¬ 
pens,  there  are  some  noticeable  differences  between  the  two 
systems.  For  the  dielectric  case,  there  is  resonant  interaction 
starting  with  1  =  4.  Also,  the  dielectric  liner  reduces  the  im- 


FIG  10  1  lie  linear  growih  rate  curves  for  a  high-energy  beam  in  a  dielec¬ 
tric  lincil  waveguide  lal  the  first  radial  mode,  ihl  the  second  radial  mode 
Parameters  n,  -  I0'‘m  ‘.ft,,  =0935,  ftL  -  0  305.  r„  =  0  06  m.  r,  =00 
m,  r4  -  0  065  m.  r„.  =  0  075  in.  and  e,  =23 
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FIG.  1 1 .  The  dispersion  curves  for  a  high-energy  beam  in  a  dielectric-lined 
waveguide:  (a)  /  =  4,  (b)  1=  15.  Parameters:  0„  =0.935,  0^  =0.305, 
r0  =  0.06  m,  r,  =  0.0  m,  r4  =  0.065  m,  rm  =  0.075  m,  and  f,  —  2.3. 


tial  harmonic  for  the  second  radial  mode  from  /  =  26  to 
1  =  9.  Furthermore,  there  is  a  dip  in  the  forward  wave 
growth  rate  at  /  =  5  for  the  lowest  curve  and  /  =  1 1  for  the 
second  curve.  These  dips  correspond  to  the  places  where  the 
resonant  points  "cross  over”  the  vacuum  light  line.  For  the 
dielectric-free  case,  the  growth  rates  for  the  forward  and 
backward  waves  are  comparable.  However,  in  the  dielectric 
case,  the  forward  wave  growth  rates  are  often  much  smaller 
than  the  corresponding  backward  wave,  especially  when  the 
resonant  intersection  is  in  region  B. 

The  dispersion  curves  for  the  fourth  and  twentieth  har¬ 
monics  are  plotted  in  Fig.  1 1;  the  beam  lines  arc  also  plotted. 
The  dispersion  curves  that  produce  resonant  intersections 
appear  to  be  approximately  hyperbolas  asymptotic  to  the 
dielectric  light  line. The  fields  are  plotted  in  Fig.  12  for  the 
backward  wave  intersection  when  /  =  4  and  for  the  second 
forward  wave  intersection  when  /  =  20.  It  is  typical  that  the 
field  is  relatively  large  at  the  beam  location  when  the  inter¬ 
section  occurs  in  region  A.  Likewise,  when  the  intersection  is 
in  region  B,  the  bulk  of  the  field  is  in  the  dielectric.  If  one 


desires  to  operate  at  a  point  in  region  B,  the  beam  must  be 
placed  very  close  to  the  dielectric  wall. 

The  effect  of  the  dielectric  is  more  pronounced  in  our 
second  example.  For  this  system,  we  have  a  low-energy  beam 
(y„  =  1  033)  with  ft,o  =  0.23,  r,  /rw  =  0.2,  and  =  0.3. 
By  plugging  these  numbers  into  Eq.  (19),  we  see  that  there 
can  be  no  resonant  interaction  for  any  /  unless  r^/r^  <0.231 
(without  some  slow  wave  structure).  Thus  for  this  example 
there  is  no  resonant  growth  at  all  without  a  dielectric  liner. 

When  we  add  a  dielectric  liner  with  a  thickness  of  3.0  cm 
and  e,  =  2.3,  we  get  resonant  interaction  starting  with 
/  =  1 1.  The  resonant  growth  rate  curve  is  plotted  in  Fig.  13. 
The  growth  rate  decreases  much  more  rapidly  with  increas¬ 
ing  /  than  the  previous  case  did.  Whereas  we  might  expect 
the  first  example  to  multimode,  the  low-energy  example  is 
likely  to  operate  only  at  the  lowest  harmonic.  The  dispersion 
curves  and  beam  line  for  the  /  =  11  case  are  plotted  in  Fig. 
14.  We  see  that  the  only  two  intersections  are  both  forward 
waves  and  the  system  is  relatively  close  to  a  tangential  inter¬ 
section  point.  Again  the  dispersion  curves  are  essentially  hy¬ 
perbolas  asymptotic  to  the  dielectric  light  line. 


Radial  Location  (m) 


FIG.  12.  The  empty  waveguide  E,  profile  for:  (a)  the  /  =  4  backward  wave 
interaction  (/-  2.94  GHz.  4,  =  2.98  m|,  (b)  the  /  =  15  forward  wav*  in¬ 
teraction  with  the  second  radial  mode  |/=  17.8  GHz,  A,  -  0014  m|.  Pa¬ 
rameters:/^  =  0.935 .0^  =  0  305.  r„  =  006  m.  r,  =  00 m,  r4  =  0065  ra. 
r_  =  0075  m,  and  r,  =  2.3. 


2876  Ptiys.FXads.VoI  28.  No.  9.  September  1985 


W.  Lawson  and  C.D. 


2876 


FIG.  1 3.  The  linear  growth  rate  curves  for  a  low-energy  beam  in  a  dielectric- 
lined  coaxial  waveguide.  Parameters:  n,  =  10u  m"J,  0a  =  0.23, 
=0.10,  r0  =  0.015  m,  r,  =0.01  m,  rd  =0.02  m,  rm  =0.05  m,  and 
f,  =  2.3. 


V.  SUMMARY 

We  have  developed  a  powerful  procedure  for  calculat¬ 
ing  the  linear  growth  rates  of  a  particular  class  of  beam- 
waveguide  systems.  The  method  easily  recovers  previous  re¬ 
sults  and  can  also  handle  systems  with  hybrid  empty 
waveguide  modes.  Using  this  method  we  have  derived  many 
properties  of  waveguide  systems  with  dielectric  liners.  We 
have  demonstrated  the  usefulness  of  such  liners  in  low-ener¬ 
gy  systems. 

Some  of  the  effects  neglected  by  our  idealized,  zero 
thickness  beam  analysis  should  be  pointed  out.  In  high  har¬ 
monic  systems  with  finite  thickness  beams,  betatron  oscilla¬ 
tions  and  frequency  (fi0)  spread  may  reduce  the  growth  rates 
considerably.  Also,  dielectric  losses  were  neglected  in  our 
analysis.  The  stabilization  of  the  negative  mass  instability  by 
lossy  dielectrics  has  been  examined  previously  by  Briggs  and 
Neil.24  These  nonideal  effects  should  be  carefully  estimated 
before  any  particular  dielectric,  large  orbit  device  is  con¬ 
structed. 

This  general  procedure  can  be  applied  to  many  more 
systems  than  the  ones  described  in  this  paper.  Although  the 
model  assumed  azimuthal  symmetry,  the  method  can  actu¬ 
ally  accommodate  systems  with  periodic  boundary  condi¬ 
tions  in  the  0  direction.  For  example,  the  method  can  be 
applied  to  the  popular  vane-type  resonator  system.25  This 
type  of  application  will  be  discussed  in  detail  in  a  future 
paper. 
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Abstract— The  production  of  high-power  microwave  radiation  from  a 
large-orbit  gyrotron  in  aaimuthally  periodic  boundary  systems  is  stud¬ 
ied  theoretically  and  experimentally.  Linear  growth  rates  are  calcu¬ 
lated  for  the  2  it  modes  of  magnetron-like  vane  resonator  (VR)  and  hole- 
and-slol  resonator  (HASR)  systems  using  a  general  growth-rate  for¬ 
malism.  The  experiment  involves  the  interaction  of  a  2.3-MeV,  1-2-kA, 
5-ns  rotating  electron  layer  with  two  ditferent  periodic  structures. 
About  500  MW  is  produced  in  Ku  band  with  a  20-slot  VR  system  and 
about  J00  MW  is  produced  in  X  band  with  a  10-slot  HASR  system.  The 
relative  merits  of  both  types  of  systems  are  discussed. 

I.  Introduction 
A.  General  Discussion 

N  RECENT  YEARS,  there  has  been  considerable  the¬ 
oretical  and  experimental  work  performed  on  the  sub¬ 
ject  of  microwave  generation  at  high  cyclotron  harmonics 
from  large-orbit  gyrotrons.  These  devices  would  be  ex¬ 
tremely  desirable  in  high-power  high-frequency  applica¬ 
tions  because  of  the  reduced  magnetic-field  requirements 
inherent  in  high-harmonic  operation.  Recent  theoretical 
studies  [1],  [2]  with  smooth  wall  geometries  indicate  that 
multimoding  would  occur  in  high-energy  systems  because 
of  the  slow  variation  of  growth  rate  with  harmonic  number 
/.  Subsequent  theoretical  studies  [3]-[7]  considered  azi- 
muthally  periodic  structures  in  attempts  to  favorably  se¬ 
lect  a  particular  harmonic.  These  studies  concentrated  on 
the  magnetron-like  vane  resonator  (VR)  structure. 

Microwave  radiation  at  cyclotron  harmonics  from  a  ro¬ 
tating  E  layer  has  been  observed  [8|  (and  subsequently- 
analyzed  [9])  in  the  Astron  machine.  More  recently,  this 
radiation  process  has  been  the  subject  of  several  experi¬ 
mental  studies  on  the  University  of  Maryland's  high-en¬ 
ergy  Rotating  Beam  Facility  (RBF)  [10|.  [II).  With  a 
smooth  cylindrical  conducting  boundary,  broad-band  ra¬ 
diation  was  observed  in  .Yband  (8-12  GHz)  at  power  levels 
of  —200  kW  per  mode.  With  a  12-vane  system  and  a  cy¬ 
clotron  frequency  of  -0.77  GHz,  single-mode  operation 
at  -  9.6  GHz  was  achieved  at  power  levels  more  than  three 
orders  of  magnitude  above  the  smooth  wall  results. 
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wj>  supported  in  pan  hy  the  Air  force  Olluc  nt  Svicnhhe  Research  and  in 
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In  this  study,  a  general  growth-rate  formalism  is  applied 
to  the  VR  configuration  and  to  the  magnetron-like  hole- 
and-slot  resonator  (HASR)  configuration.  For  both  cases, 
the  properties  of  the  dispersion  and  growth-rate  curves  are 
considered  in  detail.  Designs  for  a  20-slot  VR  structure 
and  a  10-slot  HASR  structure  are  analyzed  according  to 
the  linear  theory.  Experimental  results  for  the  two  designs 
are  then  presented  and  discussed. 

The  theoretical  formalism  is  discussed  in  Sections  I-B 
and  I-C  and  the  experimental  apparatus  is  described  in 
Section  I-D.  In  Section  II,  the  VR  structure  is  presented. 
The  HASR  structure  is  the  subject  of  Section  III  and  con¬ 
clusions  are  drawn  in  Section  IV. 

B.  Theoretical  Formalism  in  Axisymmetric  Systems 

In  this  paper,  linear  growth  rates  are  calculated  for  the 
azimuthally  periodic  structures  via  a  general  growth-rate 
formula  that  was  derived  previously  by  the  authors  [12). 
The  general  formalism  is  valid  fora  special  class  of  cylin¬ 
drical  waveguide  systems.  The  systems  have  thin,  ten¬ 
uous,  large-orbit,  annular  electron  beams.  Furthermore, 
the  systems  are  axisymmetric  and  invariant  under  trans¬ 
lations  in  the  ;  direction.  For  systems  that  satisfy  the  above 
assumptions,  an  explicit  formula  gives  the  linear  growth 
rates  of  the  EM  waves.  In  Section  1-C,  it  is  shown  that  the 
axisymmetric  growth-rate  formula  can  be  applied  directly 
to  azimuthally  periodic  systems  in  an  appropriate  limit. 

The  type  of  system  considered  hy  the  general  formalism 
has  an  electron  beam  propagating  in  some  cylindrical 
waveguide  immersed  in  a  uniform  magnetic  field  Z?.,  (see 
Fig.  1).  The  beam  is  assumed  to  be  cold:  all  electrons 
move  predominately  on  helical  orbits  about  the  system 
axis.  The  beam  is  also  assumed  to  be  very  thin;  its  equi¬ 
librium  density  is  given  by 

n0(r)  -  n Mr  -  r0)  (1) 

where  r„  is  the  equilibrium  beam  radius  and  n,  is  the  sur¬ 
face  particle  density.  Furthermore,  the  beam  is  assumed 
to  be  sutliciently  tenuous  so  that  the  dc  sell -fields  arc  ne¬ 
glected.  Each  electron  has  an  initially  unperturbed  (nor¬ 
malized)  velocity  ol  (0.  (iHi,  tf.j.  This  simple  distribution 
enables  the  use  ol  the  single-particle  equations  in  the  anal¬ 
ysis. 
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Fig.  I.  The  equilibrium  model. 

The  linear  growth-rate  problem  is  divided  into  two  parts. 
The  first  part  involves  the  calculation  of  the  perturbed 
sources  in  terms  of  the  EM  fields  via  the  equations  of  mo¬ 
tion  and  charge  conservation.  The  linearized  and  Fourier- 
analyzed  (exp  [i(kzz  +  10  -  wf)]  dependence)  equations 
are  solved  for  the  perturbed  current  to  lowest  order  in  the 
parameter  p,  =  u  -  IQ0  -  ckz(3^: 


<0^  \7o  / 


+  n0(r3rcB;)°]  S(r  -  r0)  (2b) 

illlL.  =  _ 

foC  \7o/  r0 

+  n0(r3rcB;)°l  5 (r  -  r0)  (2c) 

where  70m0c2  is  the  equilibrium  particle  energy,  fl0  =  ^B0I 
(m070),  (1.  =  kzc  -  u (3Z0,  v  =  nse2 r0n0/(2w0)  is  Budker’s 
parameter,  and  the  superscript  zeros  indicate  that  the  fields 
are  to  be  evaluated  at  the  equilibrium  radius  r0. 

The  parameter  p,  is  a  measure  of  the  difference  between 
the  frequencies  of  the  actual  EM  wave  and  the  beam  per¬ 
turbation  («0  =  /«o  +  ckzP :o).  The  synchronous  limit  is 
defined  by  \p,  -*  0.  Thus  a  beam  is  said  to  be  synchronous 
with  an  EM  wave  when  \pt  is  small  compared  to  the  other 
frequencies  in  the  analysis  (pt  «  fl0  in  particular).  A 
beam  is  said  to  be  resonant  if  is  equal  to  the  frequency 
of  an  empty  waveguide  mode.  This  paper  considers  exclu¬ 
sively  the  synchronous  resonant  case. 

In  the  second  part  of  the  analysis,  the  r  direction  is  di¬ 
vided  into  two  regions:  Region  I,  where  r  <  r0,  and  Re¬ 
gion  II,  where  r  >  r0.  In  those  regions,  7  =  0  and  the 
axial  fields  arc  written  as  a  linear  combination  of  Bessel 
functions  (the  homogeneous  wave  equation  solutions): 

X(r)  =  AJ^r)  +  BY,(Zr)  (3) 

for  =  (u/c)2  -  k2;  and  x  =  {£-.  cfl.}.  The  boundary 
conditions  that  do  not  involve  the  beam  arc  applied  in  or¬ 
der  to  find  the  solutions  in  the  two  regions  in  terms  of  the 
unknowns  £'?  and  (rdrcB.)".  The  system  of  equations  is 
closed  by  integrating  two  components  of  Ampere’s  law 
across  the  beam.  Two  equations  arc  produced  that  express 
the  jumps  in  cB.  and  rd,£.  across  the  beam  in  terms  of 


certain  beam  parameters  and  the  values  of  ( rdrcB.)°  and 

(£f. 

Although  the  beam  generally  couples  the  axial  fields  to¬ 
gether,  the  synchronous  resonant  growth  rates  for  homo¬ 
geneous  TE  modes  depend  only  on  cB.  and  its  derivative. 
The  empty  waveguide  solutions  and  corresponding  growth 
rates  for  homogeneous  TE  modes  are  given  by 


s  cB.(ro)  -  cB.(r0) 
8bb  ~  ( rdrcBf 


and 


rTF.  _  2i>lyn  U3 

Go  2  G 


(4b) 


for  T  =  Im  (w).  Equation  (4)  is  applied  to  VR  systems  in 
Section  II  and  to  HASR  systems  in  Section  III. 


C.  Extension  of  the  Growth-Rate  Formula  to 
Azimuthally  Periodic  Structures 

The  systems  considered  in  this  paper  are  assumed  to  be 
invariant  under  translations  in  the  z  direction.  Conse¬ 
quently,  the  general  EM  wave  solution  can  be  written  as 
a  sum  over  the  azimuthal  harmonic  numbers  /.  In  the  axi- 
symmetric  analysis,  all  the  boundary  conditions  are  sat¬ 
isfied  with  only  one  harmonic  perturbation  and  the  ampli¬ 
tudes  of  all  the  other  harmonic  perturbations  are  set  to 
zero.  Without  azimuthal  symmetry,  the  harmonic  pertur¬ 
bations  are  coupled  together.  Nonetheless.  (2)  is  valid  for 
each  l,  provided  that  £ ®  and  ( rdrcBz)°  are  interpreted  to 
be  the  coefficients  corresponding  to  the  / th  harmonic  of 
the  fields. 

The  synchronous  condition  implies  that  the  actual  EM 
wave  frequency  is  near  one  of  the  beam  frequencies.  Be¬ 
cause  each  beam  mode  is  separated  by  G0.  only  one  of  the 
^,’s  will  be  small  (if  p,  «  G0,  then  +  ~  «G0).  If  /o 

corresponds  to  the  synchronous  beam  mode,  then  as 
\p,a  -*  0,  7,o  becomes  much  larger  than  all  the  other  per¬ 
turbed  current  harmonics.  Consequently,  the  synchronous 
assumption  implies  that  only  one  harmonic  in  the  beam  is 
excited.  Thus  (rdr  £'.),,  and  (cB.)/n  will  be  subjected  to  jump 
conditions  at  r0,  but  all  the  other  field  harmonics  will  be 
continuous  across  the  beam.  This  implies  that  there  are 
still  only  two  boundary  conditions  that  involve  the  beam. 
Hence,  the  boundary  conditions  away  from  the  beam  can 
be  used  to  find  the  field  profiles  in  terms  of  two  unknowns. 
In  other  words,  the  linear  growth-rate  formula  can  be  ap¬ 
plied  directly  to  the  periodic  case  without  any  additional 
assumptions. 

The  solution  is  still  found  by  dividing  the  space  in  the 
waveguide  into  appropriate  regions  and  applying  all  of  the 
boundary  equations  away  from  the  beam.  One  of  the 
boundaries  is  at  the  inner  radius  of  the  slotted  wall.  The 
exact  solution  involves  matching  the  EM  fields  just  above 
and  below  this  boundary  at  every  azimuthal  location.  The 
first  simplifying  assumption  is  that  the  slots  only  contain 
the  /  =  0  harmonic.  This  means  that  the  EM  tields  remain 
constant  across  a  given  slot  opening  and  they  cannot  all 
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Fig.  2.  The  general  experimental  configuration. 


be  matched  exactly.  Instead,  one  of  the  fields  ( cB .)  is 
matched  on  the  average  across  the  opening.  This  tech¬ 
nique  corresponds  to  the  impedance  matching  technique 
used  in  early  magnetron  research  [131.  This  assumption 
enables  the  empty  waveguide  dispersion  relation  to  be 
written  in  a  reasonably  simple  form.  The  dispersion  rela¬ 
tion  can  in  principle  be  found  without  this  assumption,  but 
the  exact  solution  requires  the  calculation  of  the  determi¬ 
nant  of  an  infinitely  large  matrix.  It  has  been  shown  that 
the  approximate  results  are  close  to  the  exact  results  in 
cases  similar  to  those  considered  in  this  paper  [14). 

In  general,  there  can  be  a  uniform  phase  shift  between 
each  pair  of  neighboring  slots,  provided  that  the  total  phase 
shift  around  the  structure  is  a  multiple  of  2tr.  This  analysis 
considers  only  the  2ir  mode  (i.e.,  all  resonators  arc  in 
phase),  because  we  believe  that  it  is  the  predominant  mode 
in  our  geometry.  Our  belief  is  based  on  an  initial  condition 
calculation  [15]  which  shows  that  the  2ir  mode  is  favored 
in  systems  where  the  beam  is  symmetrically  injected  into 
the  resonator.  For  the  2ir  mode,  it  can  be  shown  that  the 
EM  fields  only  contain  multiples  of  «,  where  ft  is  the  num¬ 
ber  of  resonators  in  the  system.  The  synchronous  har¬ 
monic  is  assumed  to  be  l0  =  n. 

Finally,  it  turns  out  that  the  empty  waveguide  modes  for 
the  slotted  wall  systems  still  decouple  into  TE  and  TM 
modes.  Only  resonant  interactions  with  magnetron-TE 
modes  are  considered  in  this  paper. 

D.  Experimental  Apparatus 

A  brief  description  of  the  experimental  apparatus  is  of¬ 
fered  here;  a  more  detailed  description  can  be  found  in 
[101.  The  RBF  is  shown  schematically  in  Fig.  2.  The  beam 
is  produced  at  the  knife-edged  cathode  of  a  field-emission 
diode.  The  rotating  E  layer  is  generated  by  passing  the 
linearly  streaming  annular  beam  through  a  balanced  mag¬ 
netic  cusp  (16).  The  cusp  is  maintained  by  two  sets  of 
pancake  coils  on  opposite  sides  of  an  iron  plate.  The 
downstream  magnetic  field  is  essentially  uniform  in  the 
region  0.15  to  2  1.0  m  from  the  cusp  center.  In  this  uni¬ 
form  region,  magnetic-field  measurements  are  accurate  to 
within  I  percent.  Nominally,  the  electron  beam  has  a  0.06- 
m  radius,  an  energy  of  -2.3  MeV.  a  current  of  1-2  kA. 
and  a  post-cusp  pulse  duration  of  -5  ns. 

The  signal  is  generated  by  the  interaction  of  the  rotating 
beam  with  the  magnetron-like  structure  and  travels 
through  the  remainder  of  the  downstream  region  (total 


length  of  -3  m)  in  a  smooth  cylindrical  waveguide.  A 
cylindrical  horn  at  the  end  of  the  drift  tube  provides  the 
transition  to  free  space.  The  transition  from  the  smooth 
wall  radius  to  the  final  horn  radius  of  -0.1S  m  occurs  in 
-0.50  m.  Measurements  of  the  radiated  power  at  this 
point  indicate  that  it  is  unpolarized  and  almost  uniformly 
distributed  over  the  solid  angle  subtended  by  the  output 
horn.  The  absence  of  any  mode  pattern  is  attributed  par¬ 
tially  to  the  resolution  limit  of  the  receiving  horn  and  par¬ 
tially  to  mode  conversion  both  at  the  periodic  wall-smooth 
wall  boundary  and  in  the  output  horn.  An  open-ended 
waveguide  (WR90)  placed  >5  free-space  wavelengths 
from  the  transmitting  horn  is  used  as  the  receiving  horn 
in  the  X-band  measurements.  A  rectangular  horn  with  a 
2.3  X  10~3  m2  cross-sectional  area  placed  >10  free-space 
wavelengths  from  the  transmitting  horn  is  used  in  the  Ku- 
band  measurements. 

The  portion  of  the  radiated  power  picked  up  by  the  re¬ 
ceiving  antenna  is  transmitted  down  a  long  rectangular 
waveguide  (  —  35  m).  The  waveguide  provides  both  signal 
dispersion  and  attenuation.  At  the  end  of  this  dispersive 
line,  the  signal  is  further  attenuated  by  a  directional  cou¬ 
pler  and  a  variable  attenuator.  The  final  power  level  is 
measured  by  a  calibrated  crystal  detector  which  is  con¬ 
nected  to  an  oscilloscope  through  a  50-12  load. 

The  total  output  power  is  estimated  by  calculating  the 
fraction  of  power  intercepted  by  the  receiving  horn  and 
correcting  for  the  attenuation  suffered  by  the  signal  along 
its  path  to  the  detector  (101.  As  al1  high-power  micro- 
wave  measurements  of  this  type,  total  peak  power  is  dif¬ 
ficult  to  measure  to  an  accuracy  of  better  than  3  dB.  The 
spectrum  is  determined  by  measuring  the  arrival  time  of 
the  signal  and  using  standard  dispersive  line  theory.  Un¬ 
certainties  in  the  frequency  and  power  measurements  due 
to  the  uncertainty  in  the  measurement  of  the  arrival  time 
Can  also  be  easily  found  from  dispersive  line  theory. 


II.  The  Vane-Resonator  Configuration 


A.  Theory 

The  VR  system  is  shown  in  Fig.  3.  The  beam  radius  is 
r0,  the  inner  wall  radius  of  the  magnetron  structure  is  r„, 
the  outer  wall  radius  is  r„  and  the  number  of  slots  is  n. 
The  angle  of  the  slot  opening  is  A0  and  the  angle  of  one 
period  is  0r.  The  region  defined  by  r  <  r„  is  denoted  by 
the  interaction  region.  The  region  defined  by  r„.  <  r  <  r, 
is  called  the  slot  region. 

Using  the  fact  that  cB.  is  bounded  as  r  —  0  and  that 
only  the  nth  azimuthal  component  of  r B.  is  discontinuous 
at  r„,  it  can  be  seen  that 


cB.  =  2  (rd,eB.)) 

/a  -  OD 

+  (rd'CW 


J,(\) 


x«J;  (*J 
Jjx) 


•  10 


• nit 


(5a) 


when  0  <  r  <  r„.  and 
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cB.  =  Z  (rd,cB.)[  e‘l 


+  ( rdrcB.yn 


Qn(*-  -t0) 

Sn  (.tw,  -Tq) 


(i rdrcBfi  [4^-1  ein9 
'  UoS„(.t0,  xj_ 


when  r0  <  r  <  rw,  and 


o,epS 


ao  <  ep  <  er 


when  r„  <  r  <  rs,  where  9p  =  9  -  pOr,  0  <  p  <  n  — 
1,  $  =  [(w/c)2  -  Jkfl1'2,  x  =  r$,  etc.,  and  the  superscript 
±  refers  to  r£.  The  standard  notation  for  the  cross  prod¬ 
ucts  of  Bessel  functions  has  been  adopted  here  [  17] : 

P„(a,  b)  =  J„(a)  Yn(b)  -  Yn(a)  Jn(b)  (6a) 

Q„(a,  b)  =  Jn{a)  rn(b)  -  Yn(a)  J'n(b)  (6b) 

Rn(a.  b)  =  J'n(a)  Yn(b)  -  r„(a)  J„{b)  (6c) 

Sn(a ,  b)  =  J’(a)  K(b)  -  Y'n(a)  J'„(b).  (6d) 

The  continuity  of  £9  at  r„  requires  that 

(  A 0/0„  /  =  0 


Fig.  3.  Cross  section  of  the  20-slot  VR  structure. 

Equation  (4)  implies  that  the  magnetron-TE  empty  wave¬ 
guide  modes  are  given  by  £>(£)  =  0  and  that  the  resonant 
growth  rate  for  these  modes  is  given  by 


2»/7q 

4A 


A 9\  ^sin  in  AO  12) 
Tr)  [  nAO/2 


pCp  J'n  fa) 
L  J'ni-K)  . 


(rdrc  B.),~ 
(rd,cB.)  * 


sin  (/AO/2)  e~ 


l/n  an  integer 


Matching  the  average  value  of  cB.  across  the  slot  opening 
results  in  a  relationship  between  ( rdrcB .)*  and  (rd,cB.)\ ]. 
The  geometry  factor  ghb  can  then  be  found  from  (4)  to  be 


The  results  for  the  non  2ir  modes  arc  quite  similar.  If 
the  phase  shift  between  two  adjacent  resonators  is  as¬ 
sumed  to  be  exp  [i7„0r]  for  0  <  /„  <  n,  then  the  correct 
empty  waveguide  dispersion  relation  is  found  by  substi¬ 
tuting  /0  +  qn  for  qn  wherever  it  appears  in  (9b).  If  the 
synchronous  harmonic  is  l0.  then  the  correct  growth  rate 
is  found  by  substituting  /„  for  n  everywhere  in  (10).  Non 
2ir  mode  growth  rates  for  VR  systems  have  been  com¬ 
puted  elsewhere  (5].  Those  results  typically  indicate  that 
the  growth  rates  are  comparable  to  the  growth  rates  of  2tr 
modes. 


Xi,J„(.r,i)  S,,(t„,  .r„)  _  ( A9\  ,r„ 

0(0  77. — :  +  \~TI~ 

L  Jn(xJ  J  \VrJ 


ro  7,’ (in)  *  sin  (nAO/2) 
r„  J',MJ  .  nAO/2 


where  the  slot  term  is  given  by 


0,(0  = 


(?o(r„.  O 

50(T,lt  T,) 


and  the  interaction  region  term  is  given  by 


D,(0  = 


■  /  (  1 1,  ) 

J', 


sin  ((/itAO/2) 
qnSOI2 


0(0  =  0,(0  -  0,(0 


B.  Design 

(9a)  The  parameters  for  the  experimental  VR  structure  are 
n  —  20.  r„  =  0.0052  in.  r,  =  0.070  m.  and  SO/0,.  =  0  535. 
Plots  of  0,(0  and  0,(0  for  these  parameters  arc  shown 
in  Fig.  4.  Solutions  of  the  empty  waveguide  equation  occur 
rot  whenever  the  D,  and  D,  curves  intersect  and  are  indicated 
in  the  ligure  by  circles.  The  radial  mode  with  which  we 
expect  predominant  beam  interaction  is  also  indicated  in 
the  ligure.  The  poles  of  l),  arc  given  by  the  zeros  of /,’„(. r„) 
for  q  =  {0,  !.•••},  and  the  poles  of  0,  arc  given 
(9c)  by  the  zeros  of  S’,,  Or*.,  .r,).  When  r.  -  r_  «  r,,  D,  varies 
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Fig.  4.  The  contributions  of  the  slot  and  interaction  regions  to  the  disper¬ 
sion  function  of  the  VR  system.  The  circles  indicate  empty  waveguide 
solutions.  Parameters:  n  =  20.  rw=  0.0652  m.  r,=  0.070  m,  ±016,  = 
0.535. 


Fig  5  The  dependence  of  the  cutotf  frequency  on  the  slot  depth.  Param¬ 
eters:  n  =  20.  r„  =  0.0652,  and  ±0/0,  =  0.535. 

rapidly  compared  to  Ds  because  the  poles  of  D ,  are  far 
more  frequent  than  the  poles  of  D ,. 

Because  D ,  is  strictly  decreasing  and  D,  is  strictly  in¬ 
creasing  (on  each  domain  where  the  functions  are  defined) 

\D'(o\  =  |d;«)|  +  |o;(0|. 

We  do  not  want  to  be  near  any  of  the  poles  of  D ,  because 
T  oc  D'($)~''3.  This  generally  means  that  wc  must  have 
|D,  |  small  at  the  desired  operating  point. 

Fig.  4  can  be  used  to  determine  the  effect  of  adjust¬ 
ments  in  the  parameters  on  the  cutoff  frequencies  of  the 
various  modes  and  the  values  of  D’($)  at  those  points.  For 
example,  increasing  r,  does  not  affect  D,,  but  will  cause 
£),  to  decrease  more  rapidly,  resulting  in  lower  cutoff  fre¬ 
quencies.  Furthermore,  the  mode  trapped  between  16  and 
16.5  GHz  will  have  a  low  growth  rate  lor  all  values  of  r,. 

The  effect  of  slot  depth  on  the  cutotf  frequency  is  further 
demonstrated  in  Fig.  5.  The  sixth  through  the  ninth  radial 
modes  are  plotted  for  /t  =  20.  r„  =  0.0652  m.  and  \0I0, 
=  0.535.  The  cutoff  frequencies  of  these  modes  are  in  the 
regime  of  beam-waveguide  resonance.  The  smooth  wall 


case  is  located  at  rjrw  =  1.0.  The  numbers  in  parentheses 
next  to  each  radial  mode  number  nr  indicate  the  limiting 
smooth  wall  mode  as  r,  -*  rw.  For  example,  the  seventh 
magnet ron-TE  radial  mode  approaches  the  first  smooth 
wall  radial  mode  of  /  =  20  when  rs  -*  rw.  As  predicted, 
all  of  the  cutoff  frequencies  are  decreasing  functions  of  the 
slot  depth. 

The  EM  wave  in  the  slot  can  be  naively  pictured  as  a 
TEM  wave  traveling  down  a  shorted  waveguide.  The  \/4 
line  indicates  the  time  it  would  take  a  TEM  mode  to  travel 
up  and  down  the  slot  twice,  assuming  that  it  was  traveling 
at  the  speed  of  light.  The  corresponding  expression  for  the 
frequency  is/  =  c/[4(r5  -  rj).  The  X/4  label  comes  from 
the  fact  that  a  quarter-wavelength  standing  wave  is  ap¬ 
proximately  set  up  in  the  slot  at  these  frequencies.  For  a 
VR  system,  the  actual  X/4  frequencies  are  given  by  the 
first  zero  of  R i(xM.,  .ts).  The  frequency  of  the  naive  esti¬ 
mate  is  usually  close  to  (but  smaller  than)  the  actual  re¬ 
sult.  The  X/2  line  roughly  corresponds  to  the  frequencies 
where  £e  is  zero  at  the  slot  entrance.  The  actual  frequen¬ 
cies  are  giyen  by  the  first  zero  of  S0(xw,  x,)  (the  first  pole 
of  Ds).  At  these  frequencies,  the  standing  wave  in  the  slot 
is  a  half-wavelength.  A  X/2  point  always  has  poor  coupling 
to  the  slot  and  should  be  avoided.  At  higher  values  of 
rt/rw  (not  plotted),  good  coupling  can  be  achieved  at  fre¬ 
quencies  corresponding  to  the  3X/4  slot  mode. 

The  lowest  azimuthal  components  of  the  £e  field  for  the 
desired  operating  point  are  plotted  in  Fig.6.  For  the  inter¬ 
action  region,  the  /  =  0  contribution  is  plotted  in  Fig.  6(a) 
and  the  /  —  20  contribution  is  plotted  in  Fig.  6(b).  The 
total  £e  field  in  the  slot  is  plotted  in  both  figures  and  is 
normalized  to  unity  at  the  slot  entrance.  For  this  case 
{rjrw  =  1.074),  the  operating  frequency  is  below  the  X/4 
frequency  and  the  standing  wavelength  is  slightly  less  than 
X/4. 

'  The  bulk  of  the  /  =  20  field  component  is  near  rH  .  Be¬ 
cause  resonant  interaction  occurs  at  /  =  20.  it  is  necessary 
to  place  the  beam  near  that  wall.  Furthermore,  a  study 
considering  the  initial  condition  problem  on  the  RBF  [15) 
concluded  that  the  beam  should  be  placed  at  a  maximum 
of  the  /  =  0  component  of  Our  0.06-cm  beam  radius 
well  satisfies  both  criteria.  In  fact,  the  best  experimental 
results  for  the  VR  system  occur  when  the  downstream 
magnetic  field  is  made  a  few  percent  larger  than  the  diode 
magnetic  field.  The  effect  of  this  unbalanced  cusp  is  to 
decrease  the  average  beam  radius  [18|,  placing  the  beam 
in  an  even  more  favorable  location. 

The  magnetron-TE  dispersion  relation  in  u  —  ck.  space 
is  a  hyperbola  given  by  w;  =  w'ut  +  c'k:,  where  wlVJ,  is  the 
cutoff  frequency  of  the  mode  under  consideration  (sec  Fig. 
5  for  wllJt).  Resonances  occur  at  axial  wavenumbers  where 
ca,  (k.)  =  <j}n(k.).  In  general  there  arc  two  values  of  k.  for 
resonances;  these  are  denoted  as  a  low-frequency  reso¬ 
nance  (sometimes  a  backward  traveling  wave)  anil  a  high 
frequency  resonance  (always  a  forward  traveling  wave). 
When  these  two  roots  converge  into  one.  this  resonance  is 
denoted  as  the  grazing  condition.  A  given  beam  mode  in¬ 
teracts  with  any  radial  mode  that  has 
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Fig.  6.  The  VR  empiy  waveguide  field  ai  the  eenler  of  the  slot  opening: 
(a)  the  /  =  0  component  and  (b)  the  /  =  n  component.  The  20  slot  VR 
parameters:  n  -  20,  n,  =  7,  r„  =  0.0652.  r,  =  0.070  m,  and  a0/0r  =• 
0.535.  See  (7). 


In  high-energy  machines,  the  resonant  frequencies  arc 
usually  very  sensitive  to  small  changes  in  the  beam  energy 
and  the  applied  magnetic  field.  This  is  especially  true  for 
the  highest  radial  mode  that  satisfies  (II).  Unfortunately, 
the  highest  resonant  radial  mode  is  usually  the  desired 
mode  for  operation.  The  design  procedures  are  further 
complicated  by  the  uncertainty  in  the  RBl;  beam  energy. 

A  plot  of  growth  rate  versus /,  =  w,/2;r  at  the  resonant 
points  is  shown  in  Fig.  7  for  nominal  beam  parameters:  «, 
=  10IJ  m':.  0On  =  0.951.  tf,,  =  0.250.  and  r„  =  0.06.  The 
beam  energy  is  adjusted  slightly  so  that  equality  in  (11)  is 
nearly  achieved  for  the  highest  radial  mode  (the  so-called 
grazing  condition  |I2|).  The  growth  rate  for  the  desired 
mode  is  reasonably  large  at  -  15.2  GHz.  Thus  the  theo¬ 
retical  considerations  seem  favorable  for  this  configuration 
anil  microwave  radiation  from  a  single  radial  mode  ap¬ 
pears  likely. 


Perturbed  beam  Frequency  (n(l0  -f  ckMJlt  )/(2»l  ((IHz) 

Fig.  7.  The  linear  growth  rates  for  the  20-slot  VR  system.  Beam  parame- 
ters:  n,  =  I014  m2.  =  0.951.  <i.,  =  0.250.  r„  =  0  06  m.  The  VR 

parameters  are  given  in  Fig.  4.  The  growth  rate  at  -15.2  GHz  is  for  the 
grazing  condition  with  the  n.  =  7  mode.  The  other  low  (-)  and  high 
(  +  )  resonant  mode  interactions  are:  w_(nr  =  4)  ~  12.2  GHz:  w  (n,  = 
5)  —  12.6  GHz:  u_(n,  =  6)  —  13.3  GHz;  w.(n,  =  6)  —  17  GHz;  and 
u*(nr  =  5)  -  17.7  GHz. 

C.  Experiment 

Previous  experimental  studies  of  high-power  microwave 
generation  from  rotating  electron  beams  in  magnetron-like 
boundary  systems  have  included  configurations  with  a) 
outer  VR  boundaries  only  [11],  b)  coaxial  systems  with 
resonators  on  the  outer  wall,  c)  coaxial  systems  with  res¬ 
onators  on  the  inner  wall  [5],  and  d)  coaxial  systems  with 
resonators  on  both  walls  (in  phase  and  180°  out  of  phase) 
(5].  Previous  attempts  to  produce  microwaves  at  or  around 
the  twentieth  harmonic  of  the  electron  cyclotron  frequency 
have  shown  that  the  slot  depth  and  the  radial  position  of 
the  beam  relative  to  the  resonators  are  of  critical  impor¬ 
tance  in  maximizing  microwave  output.  Maximum  micro- 
wave  power  observed  at  the  twentieth  harmonic  using  in¬ 
ner  wall  resonators  is  about  40  MW. 

The  best  n  =  20  VR  results  to  date  were  obtained  using 
an  outer  wall  resonator  system  with  r„  =  0.0652  m,  rs  = 
0.070  m,  and  A(?/0r  =  0.535.  The  resonators  are  placed  14 
cm  back  from  the  iron  plate  so  that  they  will  not  scrape 
off  the  beam  during  its  initial  expansion.  A  typical  result 
for  the  VR  experiment  is  shown  in  Fig.  8.  Approximately 
500  MW  of  power  is  observed  at  a  frequency  of  15.5  GHz. 
Power  levels  at  other  frequencies  arc  typically  down  by  a 
factor  of  15  dB.  For  this  case,  the  diode  magnetic  field  is 
set  to  1450  G  and  the  downstream  field  is  set  to  1500  G. 
The  average  field  of  1475  G  is  very  close  to  the  cutoff  field 
of  1525  G  (see  [16)).  Consequently,  the  axial  velocity  is 
relatively  low  and  the  beam  density  is  high.  This  result 
represents  about  a  15-percent  conversion  efficiency  of 
beam  power  to  RF  power. 

The  error  bars  in  Fig.  8  do  not  reflect  the  uncertainty 
in  the  power  level  due  to  possible  systematic  error  in  the 
general  measurement  scheme  (e  g.,  error  in  the  estimation 
of  the  f  raction  of  power  intercepted  by  the  receiving  horn). 
In  Section  ID.  we  estimated  that  error  to  be  ±3  dB.  In¬ 
stead,  the  error  bars  incorporate  only  the  uncertainties 
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Fig.  8.  Radiated  microwave  spectrum  from  the  20-slot  VR  system.  Diode 
magnetic  field:  1450  G.  Downstream  field:  1500  G.  Approsimate  beam 
parameters  are  given  in  Fig.  7.  The  error  bars  incorporate  only  the  un¬ 
certainties  in  the  measurements  of  the  arrival  time  and  the  detector  volt¬ 
age.  They  do  not  reflect  possible  systematic  errors  in  the  power  mea¬ 
surement  scheme. 


in  the  measurements  of  the  arrival  time  and  the  detector 
voltage.  In  this  way,  relative  comparisons  between  the 
results  of  the  various  wall  configurations  can  reasonably 
be  made.  Unfortunately,  the  large  net  error  estimate  ( ±  3.5 
dB)  leaves  the  uncertainty  in  the  efficiency  to  be  (more 
than)  a  factor  of  two. 

The  diagnostics  for  the  experiment  leave  several  ques¬ 
tions  unanswered.  For  example,  at  these  high  harmonics, 
the  mode  spacing  is  quite  small  and  the  uncertainty  in  the 
beam  energy  makes  it  impossible  for  us  to  identify  the  ex¬ 
act  mode  from  the  frequency  measurements.  Thus  it  is 
possible  that  the  machine  is  actually  operating  in  a  mode 
near  the  synchronous  2jt  mode.  The  experimental  results 
are  consistent  with  the  27r  mode  theory,  but  they  do  not 
confirm  it. 

In  spite  of  the  uncertainties,  several  conclusions  can  be 
drawn  from  this  and  previous  experiments  on  the  RBF. 
First,  a  properly  designed  periodic  structure  can  excite 
radiation  at  a  single  frequency  with  power  levels  three  or¬ 
ders  of  magnitude  above  those  from  a  smooth  wall  config¬ 
uration.  Also,  the  output  from  this  20-slot  VR  experiment 
is  comparable  in  power  to  the  best  12-slot  VR  experiments 
previously  performed  on  the  RBF.  This  result  indicates 
that  the  outlook  for  operating  at  extremely  high  cyclotron 
harmonics  is  promising  (although  initial  30-slot  VR  ex¬ 
periments  had  severe  multimoding  problems  and  achieved 
power  levels  at  best  a  factor  of  10  higher  than  smooth  wail 
results  (I9|).  Finally,  it  should  be  pointed  out  that  in  pre¬ 
vious  studies  with  periodic  structures  on  the  inner  wall 
w.  radiation  in  the  frequency  range  consistent  with  the 
ir  mode  remained  at  smooth  wall  power  levels  (were  not 
enhanced)  in  systems  that  radiated  effectively  in  the  fre¬ 
quency  range  consistent  with  the  2t  mode.  This  fact  lends 
some  support  to  the  initial  condition  calculation,  because 
the  cusp  injects  the  rotating  beam  symmetrically  into  the 
region  with  the  periodic  structures. 


Fig.  9.  Cross  section  of  the  10-slot  HASR  structure. 


III.  The  Hole  and-Slot  Configuration 
A.  Theory 

The  HASR  system  is  shown  in  Fig.  9.  It  is  assumed  that 
the  HASR  system  consists  of  a  VR  system  with  a  hole 
added  at  the  end  of  each  vane.  The  rectangular  slot  case 
is  similar  and  is  not  considered  here.  The  hole  radius  is  rh 
and  the  angle  of  the  hole  opening  is  Ai p  =  2  sin'1  [sin 
(A0/2)  rjrh). 

An  additional  assumption  has  to  be  made  for  the  HASR 
system.  Let  r  and  0  denote  the  cylindrical  coordinates 
whose  origin  is  at  the  center  axis,  and  let  r*  and  0*  denote 
the  cylindrical  coordinates  whose  origin  is  at  the  center  of 
a  (hole)  resonator.  At  the  opening  of  the  hole,  the  arc  de¬ 
fined  by  |r|  =  r,  is  assumed  to  coincide  with  the  arc  de¬ 
fined  by  |r*|  =  rh.  This  assumption  will  be  good  only  for 
relatively  thin  slots.  Mathematically,  the  restriction  is  that 
A 0  «  rh/rs.  This  assumption  simplifies  the  calculation 
of  the  boundary  conditions  at  the  hole  opening,  because 
the' fields  in  the  resonator  are  easily  expressed  in  terms  of 
the  (r*.  9*)  coordinate  system. 

The  field  in  the  slot  is  now 


cB.  =  ( rdrcB .)* 


Qi i(t.  -ri) 

jck5o(.v„.  x,) 


+  (rd,B:)' 


(?o(t.  O 
*A(ri.  xj 


(12a) 


and  the  field  in  the  hole  resonator  is 


cB, 


( rd,cB .) 


v;cu) 


c 


(12b) 


where  r*  =  r** .  .r,,  =  r;,£,  and  ihc  field  in  the  interaction 
region  is  still  given  by  (5a)  and  (5b). 

The  requirement  (hat 
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comes  from  matching  £e  and  the  average  cB.  field  at  r  = 
r,.  At  this  point,  the  solution  can  be  found  by  following 
the  steps  for  the  VR  solution,  because  there  is  only  one 
remaining  unknown  in  the  slot  and  all  the  boundary  con¬ 
ditions  at  rw  remain  to  be  applied. 

The  dispersion  relation  can  be  split  into  three  parts.  The 
slot  term  and  the  interaction  term  are  still  given  by  (9a) 
and  (9b).  The  hole  resonator  term  is  given  by 

Dh'iZ)  =  x»xsS0(xw,  xs) 


X 


At L 

*s)  +  So(-r» 
2ir 


xs) 


sin  (mA\i/2) 


mA^/2  J 


04) 


The  equations  for  the  geometry  factor  (8)  and  the  growth 
rate  (10)  are  valid  for  the  HASR  system  after  the  substi¬ 
tution 


0(0  =  0,(0  +  0,( sO  -  0,(0.  (15) 

It  is  often  more  convenient  to  plot  the  entire  resonator  term 
0,(0  -  0,(0  +  0,(0-  From  (9a)  and  (14),  it  is  easy  to 
show  that 


Fig.  10.  The  contributions  of  the  resonator  and  interaction  regions  to  the 
dispersion  function  of  the  HASR  system.  The  circles  indicate  empty 
waveguide  solutions.  Parameters:  n  =  10,  -  0.068  m.  r,  =  0.070  m. 

rh  =  0.012  m.  ±0/0,  =  0.147. 

est  radial  mode  competitor  may  not  be  as  complete  with 
VR  systems.  The  added  flexibility  in  the  choice  of  rh  is 
another  advantage  of  HASR  systems. 

In  Fig.  11,  the  total  £v  field  profile  of  the  fifth  radial 
mode  is  plotted  as  a  function  of  r.  The  field  along  the  line 


Po(.xw,  xs)  +  (A^I2ir)  Q0(xw,  xs)  £  JJx,,VJ’Jxh)  [sin  (wAv^/2)/mA^/2): 

0,(0  =  - — - .  (16) 

xs)  +  (A\J/2t)  Su(xw,  xs)  £  Jjxh)/J'm(xh)  [sin  (mA^I2)/tnA^I2]: 


B.  Design 

The  parameters  for  the  HASR  system  are  n  =  10,  rw  = 
0.068  m.  r,  =  0.0703  m,  rh  =  0.012  m,  and  A0ldr  =  0.147. 
The  dispersion  functions  Dr  and  D,  arc  plotted  in  Fig.  10. 
Recall  that  solutions  are  given  by  the  intersections  of  the 
two  curves.  In  the  figure,  these  points  arc  indicated  by 
circles.  The  system  is  designed  to  operate  in  the  fifth  ra¬ 
dial  mode.  The  cutolf  frequency  for  that  mode  is  7.63  GHz. 
The  poles  and  the  zeros  of  Dr  can  be  found  from  (16). 

The  analogy  of  TEM  X/4,  X/2,  etc.  modes  is  more  ab¬ 
stract  for  the  HASR  resonators,  but  still  can  be  applied. 
The  X/2  point,  for  example,  corresponds  to  the  first  pole 
of  Dr  (16).  In  the  experimental  design,  the  X/2  point  is  at 
~7  GHz  and  the  desired  operating  point  is  between  that 
and  the  3X/4  point.  The  proximity  of  the  X/2  point  to  the 
mode  at  7.13  GHz  results  in  a  large  value  of  |D'($)I  and 
consequently  a  small  value  of  growth  rate  (see  ( 10)).  How¬ 
ever.  D,  decreases  sufficiently  rapidly  so  that  the  desired 
mode  has  a  favorable  growth  rate.  Thus  careful  placement 
of  the  first  pole  of  Dr  can  virtually  eliminate  competition 
from  the  radial  mode  nearest  the  desired  operating  point. 

The  20-slot  VR  case  previously  discussed  was  designed 
for  the  X/4  mode,  had  slowly  varying  growth  rates,  and 
should  not  be  compared  with  tins  HASR  example.  How¬ 
ever.  the  same  idea  can  be  applied  to  the  VR  systems  op¬ 
erating  between  the  X/2  and  3X/4  points.  In  general,  l)r 
varies  more  rapidly  than  D,  and  so  the  nulling  of  the  near¬ 


passing  through  the  center  of  the  slot  is  shown  in  Fig. 
11(a).  The  £„  field  midway  between  slots  is  plotted  in  Fig. 
11(b).  The  beam  is  located  in  a  reasonably  favorable  po¬ 
sition  of  the  field  profile.  At  small  radii,  the  /  =  0  term 
dominates  and  the  field  is  essentially  independent  of  azi¬ 
muthal  angle.  However,  near  rlv,  the  /  =  n  term  is  quite 
strong  and  the  £#  field  is  essentially  n  periodic.  Compared 
to  the  VR  example,  the  relative  contribution  of  the  /  =  0 
term  is  greatly  reduced.  This  is  partly  due  to  the  smaller 
value  of  AOIOr  and  partly  due  to  the  lower  radial  mode 
number.  The  £t  field  in  the  hole  is  quite  large,  though  the 
hole  volume  is  relatively  small. 

In  the  experiment,  the  best  results  occurred  at  magnetic 
fields  between  1300  and  1325  G.  The  growth  rates  are 
shown  in  Fig.  12  for  the  nominal  parameters  dAl  =  0.854. 

=  0.487.  and  r„  =  0.06  in.  The  growth  rates  for  the 
fifth  radial  mode  (at  -8.1  and  —9.7  GHz)  arc  approxi¬ 
mately  six  times  larger  than  the  growth  rates  of  the  fourth 
radial  mode  (at  -7.2  and  -10.6  GHz),  indicating  that 
single-mode  operation  should  be  feasible. 

C.  Experiment 

The  HASR  system  is  inserted  into  the  RDF  and  held  in 
place  -15  cm  from  the  iron  plate.  An  adapter  brought  the 
wall  radius  down  from  the  RDF's  7.5  cm  to  the  HASR'x 
6.8  cm.  The  overall  length  of  the  HASR  structure  is  ap- 
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Fig  II  The  HASH  empty  waveguide  £,  field:  (a)  at  the  slot  center  and  (b) 
midway  between  two  slot  centers.  The  lO-slot  HASR  parameters,  n  = 
10.  n.  =  5 ,  r„  =  0  068  m,  r,  =  0  070  m,  r„  =  0  012  m.  and  A 6/9 ,  = 
0.147. 


Fig  12  The  linear  growth  rates  for  the  10  slot  HASR  system  Hearn  pa¬ 
rameters  n,  =  It)'*  m  =  0  854.  —  0  487.  and  r„  -  00b  m. 

The  HASR  parameters  are  given  in  Fig  10.  The  mode  identifications 
are  w.ln,  =4)  —  72(ill/.uj  In,  =  5|  -  R  1  Gllr.  u.tn,  “  5)  -  9.7 
Oil/,  and  ut.ln,  -4)  —  100(111/ 

proximatcly  60  cm.  Short  rods  with  20°  tapers  arc  in¬ 
serted  in  the  holes  at  the  input  side  ot  the  structure  to  pro¬ 
vide  a  continuous  transition  trom  the  smooth  wall  conlig- 
uration  to  the  periodic  configuration.  The  MASK  structure 


Fig.  13.  Radiated  microwave  spectrum  from  the  10-slot  HASR  system 
Magnetic  held:  1300  G.  Approximate  beam  parameters  are  given  in  Fig. 
12.  The  error  bars  incorporate  only  the  uncertainties  in  the  measurements 
of  the  arrival  time  and  the  detector  voltage.  They  do  not  reflect  possible 
systematic  errors  in  the  power  measurement  scheme. 

is  tapered  at  the  output  end  so  that  the  large  amount  of 
field  energy  in  the  hole  can  easily  be  collected.  The  24° 
taper  brings  the  inner  radius  out  to  a  smooth  wall  radius 
(rw  >  rs  +  2 rh).  A  10-cm-radius  smooth  wall  tube  brings 
the  microwave  radiation  from  the  end  of  the  HASR  system 
to  a  cylindrical  cone  at  the  end  of  the  vacuum  chamber. 

The  best  results  for  this  experiment  are  shown  in  Fig. 
13.  Approximately  260  MW  of  power  is  observed  at  a  fre¬ 
quency  of  9.2  GHz  in  a  balanced  field  of  1300  G.  Power 
levels  at  other  frequencies  are  typically  down  by  more  than 
10  dB.  Single-mode  operation  is  consistently  achieved  at 
power  levels  exceeding  200  MW.  By  comparing  Figs.  12 
and  13,  it  appears  these  conditions  represent  the  forward- 
wave  resonance  with  the  nr  =  5  waveguide  mode. 

It  is  unfortunate  that  the  sharp  resonance  at  1300  G  oc¬ 
curred  so  far  below  the  cutolf  field  (1525  G).  The  beam 
density  at  this  point  is  only  one-half  of  the  density  at  the 
20-slot  VR  resonance  point  (because  /?.„  at  1300  G  is  nearly 
twice  the  value  of  /3%l  at  1475  G).  This  reduction  in  n,  prob¬ 
ably  accounts  (at  least  partially)  lor  the  reduced  power  lev¬ 
els  in  the  10-slot  HASR  device. 

The  error  bars  on  the  frequency  measurement  result 
from  the  uncertainty  in  the  arrival-time  measurement. 
Again,  the  power-level  error  bars  reflect  only  the  uncer¬ 
tainties  in  the  arrival  time  and  the  detector  voltage. 

IV.  Summary 

In  this  paper,  wc  have  show  n  that  the  results  of  a  general 
axisymmetric  linear  growth-rate  formalism  arc  applicable 
to  azimuthally  periodic  systems  in  the  synchronous  limit. 
We  used  the  growth-rate  formula  to  calculate  the  proper¬ 
ties  of  magnetron  like  VR  and  HASR  structures.  Wo  lound 
that  the  HASR  system  theoretically  has  an  advantage  over 
the  VR  system  in  that  its  parameters  can  he  more  readily 
adjusted  to  suppress  radial  mode  competition.  Unfortu¬ 
nately.  an  additional  assumption  in  the  HASR  analysis 
forced  a  restriction  on  the  maximum  allowable  size  for  the 
relative  slot  opening  (S0/0,)-  This  is  undesirable  he- 
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cause  of  the  resultant  reduction  in  slot  coupling 

[r  OC  (A0/0,)1'3]. 

Experiments  with  a  20-slot  VR  system  and  a  10-slot 
HASR  system  successfully  produced  high-power  micro- 
waves  at  power  levels  comparable  to  previous  12-slot  VR 
systems.  The  radiation  in  each  case  was  determined  to  be 
at  a  single  frequency  (corresponding  to  a  single  mode)  to 
within  the  accuracy  of  the  dispersive  line  measurements. 
Both  systems  had  relatively  stable  operation  at  frequencies 
consistent  with  synchronous  2v  modes,  although  the  modes 
were  not  positively  identified  as  such.  The  results  indicate 
that  both  of  these  configurations  may  be  well  suited  to 
high-power  generation  in  the  centimeter  and  millimeter 
regime. 

We  have  successfully  demonstrated  that  more  than  one 
type  of  configuration  can  be  used  to  promote  single-mode 
operation  in  high-power  large-orbit  gyrotrons.  However, 
much  work  remains  to  be  done  on  these  configurations  be¬ 
fore  optimal  designs  can  be  built.  Alternate  schemes  for 
periodic  structures  can  be  tested  with  the  linear  theory  via 
(4).  Nonlinear  theories  and  particle  simulations  need  to 
address  questions  of  efficiency.  Other  experiments  should 
be  performed  to  determine  the  actual  operating  mode,  test 
the  initial  condition  calculation  [15],  and  evaluate  the  suit¬ 
ability  of  various  resonator  configurations. 
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The  generation  of  millimeter  wave  radiation  from  the  interaction  of  a  rotating  electron  beam  (2 
MeV,  1  kA,  5  nsec)  with  an  azimuthally  periodic  wiggler  magnetic  field  has  been  studied 
experimentally.  Calculations  of  the  effects  of  the  wiggler  magnetic  field  on  the  single  particle 
electron  orbits  are  presented,  together  with  experimental  measurements  of  the  effects  of  the 
wiggler  field  on  the  electron  beam.  Narrow -band  radiation  at  power  levels  in  excess  of  200  k  W  has 
been  observed  at  88  and  175  GHz  for  wiggler  fields  with  6.28  and  3. 14  cm  periods,  respectively. 

The  radiation  frequency  spectra  for  various  experimental  configurations  are  presented,  and 
results  are  compared  with  theoretical  expectations. 


I.  INTRODUCTION 

In  recent  years  many  theoretical  and  experimental 
studies  have  been  reported  of  linear  free-electron  lasers 
(FEL's)  in  which  short-wavelength  radiation  is  produced  by 
the  interaction  of  an  electron  beam  with  a  spatially  periodic 
wiggler  magnetic  field. 10  Recently,  a  novel  circular  geome¬ 
try  FEL  has  been  explored  both  theoretically  and  experi¬ 
mentally  in  a  collaborative  effort  by  researchers  at  the  Mas¬ 
sachusetts  Institute  of  Technology  and  the  University  of 
Maryland."'16  In  this  concept,  a  rotating,  relativistic  elec¬ 
tron  beam  interacts  with  an  azimuthally  periodic  wiggler 
field  produced  by  samarium  cobalt  magnets  placed  interior 
and  exterior  to  the  beam.  The  potential  advantages  of  such 
systems  include  a  longer  effective  interaction  region,  a  more 
compact  geometrical  configuration,  and  internal  feedback 
resulting  from  the  recirculation  of  the  electromagnetic  wave. 
This  last  feature  may  mean  that  the  device  can  operate  as  an 
oscillator  rather  than  an  amplifier,  as  in  the  case  of  linear 
FEL’s. 

In  the  experiments  to  date,  two  principal  methods  have 
been  used  to  generate  the  rotating  electron  beam.  The  first 
experiments  employed  a  diode  configuration  similar  to  those 
used  in  relativistic  magnetrons."14  Here  the  electrons  per¬ 
form  ExB  drifts  around  the  aximuth  in  the  presence  of  a 
radial  electric  field  and  an  axial  magnetic  field.  Addition  of 
an  aximuthally  periodic  magnetic  field  then  results  in  a  cir¬ 
cular  FEL.  Although  initial  experimental  results  from  this 
configuration  are  encouraging,  one  potential  drawback  of 
this  configuration  is  the  considerable  electron  velocity  shear 
inherent  to  cross  field  electron  flow  in  magnetron-like  de¬ 
vices. 

Recently  we  published  a  preliminary  study  of  a  second 
configuration  that  effectively  circumvents  this  velocity  shear 
problem.15  This  experiment  involves  the  generation  of  an 
essentially  monocnergetic  rotating  electron  beam  by  passing 
a  hollow,  nonrotating  electron  beam  through  a  narrow  mag¬ 
netic  cusp.  In  this  manner,  the  v,  X B,  force  at  the  center  of 
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the  cusp  effectively  converts  the  axial  beam  velocity  into  ro¬ 
tational  velocity  downstream  of  the  cusp  transition  region.  If 
the  cusp  is  symmetric,  the  downstream  beam  performs  sim¬ 
ple  axis  encircling  cyclotron  orbits  with  a  gyroradius  equal 
to  the  radius  of  the  beam  on  the  upstream  side  of  the  cusp.17 

To  the  axial  magnetic  field  about  which  the  beam  elec¬ 
trons  rotate  is  added  an  azimuthally  periodic  wiggler  field, 
which  is  primarily  radial  near  the  center  of  the  gap,  and  thus 
transverse  to  the  electron  beam  flow.  The  wiggler  field  is 
produced  by  samarium  cobalt  magnets  placed  behind  two 
concentric  metal  cylinders  in  such  a  way  that  the  beam  sees 
only  smooth  conducting  boundaries.  The  proximity  of  the 
conducting  walls  to  the  beam  also  serves  to  suppress  the 
negative  mass  instability,  which  has  been  used  to  produce 
radiation  at  microwave  frequencies  in  previous  studies.1* 19 

In  this  paper,  we  present  the  first  detailed  measure¬ 
ments  of  the  operational  characteristics  of  this  new  circular 
FEL  configuration.  Measurements  of  the  effect  of  the 
wiggler  field  on  the  rotating  electron  beam  are  presented,  as 
are  radiation  spectra  obtained  for  a  number  of  wiggler  con¬ 
figurations.  Section  II  of  this  paper  contains  a  discussion  of 
theoretical  considerations,  and  the  experiments  are  present¬ 
ed  in  Sec.  III.  Conclusions  are  drawn  in  Sec.  IV. 


II.  THEORETICAL  DISCUSSION 

A.  Electron  motion  in  the  axial  and  wiggler  magnetic 
fields 

The  general  configuration  used  for  these  studies  is  de¬ 
tailed  in  Fig.  I .  Downstream  of  the  cusp  transition,  the  elec¬ 
tron  orbits  have  axial  (ul2 )  and  azimuthal  (uft2 )  velocity  com¬ 
ponents  given  in  terms  of  the  upstream  axial  velocity  (u,, )  by 

*£i  =  +"«:  =^2  + 

where  rc  is  the  cathode  radius  and  /?,,  =  t is  the 
relativistic  electron-cyclotron  frequency  in  the  downstream 
axial  magnetic  field.  Thus,  as  the  magnetic  field  is  raised 
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that  near  the  center  of  the  gap,  the  field  is  primarily  radial 
and  that  the  undulatory  v<,  xB,  force  is  the  ±  z  direction, 
analogous  to  the  transverse  motion  of  electrons  in  a  linear 
FEL. 

Because  of  the  complicated  nature  of  the  combined  axi¬ 
al  and  wiggler  fields,  a  single  particle  computer  simulation 
program  has  been  used  to  check  that  the  particle  orbits  are  as 
desired.  Figure  2  shows  calculated  particle  orbits  with  and 
without  the  wiggler  fields,  results  that  show  clearly  that  the 
electron  orbits  are  almost  unperturbed  in  the  r-6  plane,  and 
the  ui.dulation  is  primarily  axial,  as  desired.  While  Fig.  2 
shows  typical  electron  orbits  for  an  electron  with  an  initial 
radial  position  of  precisely  6  cm  (the  cathode  radius),  elec¬ 
trons  launched  at  other  radial  positions  within  the  radial 
width  of  the  beam  do  not  show  significantly  different  behav¬ 
ior.  We  have  performed  these  calculations  for  each  of  the 
several  wiggler  field  configurations  investigated  experimen¬ 
tally,  and  in  no  case  are  the  electron  orbits  unsatisfactory. 
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B.  Excitation  of  TM  waves  by  the  rotating  electron  beam 

In  a  previous  theoretical  paper  by  two  of  the  authors  (Y. 
Yin  and  G.  Bekefi),"*  the  radiative  process  has  been  identi¬ 
fied  as  the  coupling  of  a  “synchronous  mode,"20  upshifted  in 
frequency  by  the  wiggler  periodicity  N-. 

(o  =  (I  +  A')/?,, 

to  one  or  more  of  the  TM  waves  supported  by  the  coaxial 
conducting  boundary  system.  The  radiation  frequency  is 
given  by 

Nil,,  _  Ar,„t>u 

1  -[//?., /wr(/,mj]  \-{v0/vp0) 


toward  a  cutofT  value,  given  by 
Blc  =  m0y/rce, 


the  axial  velocity  is  reduced  and  electron  orbits  become  fair¬ 
ly  tight  helices.  Beyond  the  cusp  the  electrons  move  in  the 
presence  of  the  combined  axial  and  wiggler  magnetic  fields, 
which  can  be  approximated  (subject  to  the  condition  that 
V-B  =  VxB  =  0  in  the  region  between  two  conducting 
boundaries  of  radii  r0  and  r , )  by  the  expression 


' r  \S  *  'l/r  V-V*  -  1 1/2* 

.7)  ](?) 


+  zB,. 


Here  r,  0,  and  z  are  unit  vectors  in  the  radial,  azimuthal,  and 
axial  directions,  respectively,  , V  =  n{r„  +  r,)//„  is  the  num¬ 
ber  of  spatial  periods  around  the  azimuth,  /„  is  the  linear 
periodicity  specified  midway  in  the  gap,  and  B„w  is  the  am¬ 
plitude  of  the  radial  component  of  the  field  at  a  distance 


,  _  /-.V  •  I  V  ,  I  l/U.Vl 

'  ~  V  ro  )  I 

where  the  azimuthal  component  vanishes.  It  is  easily  seen 


(b) 


FIG  2  Calculated  particle  orbits  in  the  r-0  and  r-i  planes  for  an  electron 
injected  with  u,  =  0  20c,  =0  df >r  into  the  interaction  space  with  lal 

B,  =  1.4  kU.  =  0,  and  (hi  H,  ~  1.4  kti,  /?„„  =  U  kU. 
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TABLE  I  Computed  parameter*  for  an  electron  ring  of  radius  6  000  cm 
rotating  in  a  coaxial  waveguide  of  radii  ra  -  b.5(N  cm.  r,  =  5.3^$  cm. 
0a  —  0.*>68;  =  0.75  k G.yrt  •=  HO  A  cm  u  is  the  radiation  frequen¬ 

cy;  lo,  the  temporal  growth  rate;  and  if  the  saturated  efficiency. 


N 

k 

|cm| 

m 

/ 

ut/2r 

|GHz| 

tit,  x  10  "* 
{rad 
sec  '| 

V 

rd 

6 

6.28 
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29 

26.8 

1.02 

19 

6 

6.28 
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94  7 
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3.14 
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14  6 

0  48 

14 

12 

3.14 

i 

281 

225.6 

0014 

0.15 

12 

3.14 

2 

56 

52.3 

0.15 

0.83 

12 

3.14 

2 

232 

187  9 

0  031 

0.05 

where  kw  =  2 r//0,  eoc{l,m )  is  the  cutoff  frequency  for  the 
TM,m  mode,  /  and  m  are  the  azimuthal  and  radial  wavenum¬ 
bers,  respectively,  and 

V  =  coc(l,m)(r/l) 

is  the  azimuthal  phase  velocity  of  the  rf  perturbation.  It  is 
interesting  to  compare  this  result  to  the  dispersion  relation 
for  a  conventional  linear  free-electron  laser: 

[1  -  (t’o/c)  ]  ’ 

where  is  now  the  axial  electron  velocity. 

In  the  limit  where  the  gap  between  the  two  coaxial  con¬ 
ducting  boundaries  is  small  compared  to  their  mean  radius, 
the  cutoff  frequency  of  the  TM,m  mode  may  be  approximat¬ 
ed  by  the  expression 

=  [ mrrc/r,{g -  1)](1  -fa/')'/2, 
where  g  =  r^/r,  and 


The  resultant  predicted  radiation  frequencies  for  actual  ex¬ 
perimental  parameters  discussed  in  the  next  section  are  sum¬ 
marized  in  Table  I.  For  the  specified  parameters  these  are 
the  only  unstable  interactions  predicted  by  the  theory.  Be¬ 
cause  these  predicted  radiation  frequencies  are  quite  sensi¬ 
tive  to  the  values  chosen  for  the  applied  magnetic  field  and 
the  electron  energy,  the  values  indicated  in  Table  I  must  be 
treated  as  estimates.  One  important  conclusion,  however,  is 
that  high-frequency  radiation  is  predicted  for  very  low  val¬ 
ues  of  the  radial  mode  number  m,  an  important  result  be¬ 
cause  the  electron  beam  radial  width  (about  5  mm)  is  a  siza¬ 
ble  fraction  of  the  gap  between  the  two  coaxial  conductors 
(13  mm).  A  comparison  between  these  predicted  values  and 
experimental  measurements  will  be  made  in  Sec.  IV.  The 
growth  rate  of  the  instability  o,  was  calculated  for  the  case  of 
an  azimuthal  current  density  J„  —  1.1  A  cm  2  in  order  to 
satisfy  the  assumption  of  a  tenuous  electron  ring,  and  is  also 
indicated  in  Table  I.  In  the  actual  experiments,  the  azi¬ 
muthal  current  density  is  estimated  to  be  —  bO  A  cm  2.  In 
the  so-called  "single  particle,  high-gain  strong  pump"  re¬ 
gime*  the  instability  growth  rate  is  proportional  to  J'„'\ 
which  implies  that  the  experimental  growth  rates  would  be 
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about  a  factor  of  4  higher  than  those  given  in  the  table.  The 
last  column  of  Table  I  lists  the  estimated  saturation  effi¬ 
ciency  t)  caused  by  phase  trapping  of  electrons  in  the  poten¬ 
tial  wells  of  the  ponderomotive  potential. 

III.  EXPERIMENTS 
A.  Apparatus 

The  general  experimental  configuration  is  shown  in 
Fig.  1.  A  hollow,  nonrotating  electron  beam  is  field  emitted 
from  a  6  cm  radius,  knife-edge  cathode  located  7.5  cm  up¬ 
stream  of  a  brass  anode  plate.  A  5  mm  wide  aperture  slit  in 
the  anode  plate  allows  a  fraction  of  the  electron  beam  cur¬ 
rent  to  pass  through  the  anode  plane  into  the  cusp  transition 
region.  The  cusp  magnetic  fields  are  produced  by  two  inde¬ 
pendently  controllable  solenoids,  and  a  soft  iron  plate  is  used 
to  narrow  the  axial  extent  of  the  cusp  transition  region.  The 
measured  FWHM  of  the  radial  component  of  the  magnetic 
field  at  the  center  of  the  cusp  is  4  mm.  The  v,  X  Br  force  acts 
to  convert  axial  electron  velocity  upstream  of  the  cusp  tran¬ 
sition  region  into  azimuthal  velocity  downstream,  with  re¬ 
sulting  downstream  beam  parameters  of  2  MeV,  1  kA,  and  5 
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FIG  .V  Arrangnient  ul  bar  magnets  (lop),  Mali  probe  measurement  of  the 
wigylcr  field  at  a  radial  position  r  —  5  92  cm.  as  a  function  of  azimuthal 
angle  Iboltomj 
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FIG.  4.  Peak  current  exiting  the  interaction  region  versus  applied  cusp  axial 
magnetic  field.  Results  with  and  without  the  wiggler  magnets  are  shown. 


nsec.  Because  the  cusp  reflects  all  electrons  with  energies  less 
than  a  threshold  value  given  by 

£.h  =  [("V2)2  +  (cercB)2] 1/1  ~  "V2. 
the  total  energy  spread  in  the  downstream  rotating  beam  is 
in  the  range  l%-3%.  The  velocity  spread  caused  by  local 
temperature  in  the  beam  has  been  estimated  to  be  about 
0.2%. 

The  downstream  beam  rotates  between  two  concentric 
stainless  steel  cylinders  of  radii  5.40  and  6.51  cm,  respective¬ 
ly.  Single  turn  Rogowski  coils  can  be  located  at  the  upstream 
and/or  downstream  ends  of  the  interaction  region  to  mea¬ 
sure  the  axial  current  entering  or  leaving  the  wiggler  region 
with  or  without  the  wiggler  magnets  in  place.  The  samarium 
cobalt  magnets  used  to  provide  the  wiggler  field  are  placed 
behind  the  cylinders  and  held  in  place  by  grooved  aluminum 
holders.  Typical  arrangements  of  the  magnets  to  achieve 
various  periodicities  are  shown  in  Fig.  3.  The  axial  length  of 
the  wiggler  field  is  about  20  cm.  and  the  wiggler  strength  has 
been  varied  in  the  experiments  by  simply  removing  some  of 
the  magnets  from  the  6.28  cm  wiggler  configuration  shown 
in  the  figure. 

The  radiation  generated  in  the  experiments  is  measured 
with  a  small  horn  antenna  located  immediately  downstream 
of  the  interaction  region.  The  radiated  power  spectra  for  var¬ 
ious  configurations  has  been  measured  using  a  sensitive  grat¬ 
ing  spectrometer,2 1  with  gratings  available  in  the  range  70- 
200  GHz.  The  frequency  resolution  of  the  spectrometer  is 
typically  4/7/  =  0.02,  and  the  insertion  loss  is  in  the  range  3- 
5dB. 


B.  Electron  beam  measurements 

Measurement  of  the  axial  electron  beam  current  exiting 
the  interaction  region  with  and  without  the  wiggler  magnets 
in  place  has  been  made  using  a  single  turn  Rogowski  coil 
located  immediately  downstream  of  the  wiggler  region.  The 
results,  shown  in  Fig.  4,  indicate  that  a  wiggler  field  of  1 300 
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FIG.  5.  Oscilloscope  waveforms  of  (a)  diode  voltage,  (b)  axial  current  col¬ 
lected  by  a  2.24  mm"’  collector  located  in  the  center  of  the  interaction  region, 
(c)  microwave  signal  inT  band  (9 1-170  GHz)  with  wiggler  magnets,  and  Id) 
microwave  signal  in  T  band  without  wiggler  magnets. 


G  amplitude  and  /„  =  6.28  cm  causes  a  drop  in  the  beam 
current  of  only  about  30%.  These  results  confirm  that  the 
wiggler  field  does  not  have  a  disastrous  effect  on  the  electron 
orbits,  a  result  consistent  with  the  single  particle  orbit  calcu¬ 
lations.  In  fact,  it  is  likely  that  the  wiggler  field  acts  to  re¬ 
move  electrons  with  badly  off-centered  orbits  from  the 
beam. 

Another  indication  of  the  nature  of  the  beam  current 
pulse  shape  has  been  obtained  by  placing  a  small-area  axial 
current  collector  midway  between  the  conducting  cylinders 
at  the  axial  center  of  the  wiggler  field  region.  A  typical  cur¬ 
rent  pulse  waveform  is  shown  in  Fig.  5  and  shows  dramati¬ 
cally  the  shortening  of  the  downstream  electron  beam  cur¬ 
rent  pulse  duration  caused  by  the  reflection  of  all  low-energy 
electrons  at  the  cusp  transition  region.  This  independent 
measurement  of  the  beam  current  in  the  wiggler  region  also 
shows  that  even  a  relatively  strong  wiggler  field  does  not 
disastrously  disrupt  the  rotating  electron  beam. 

C.  Radiation  measurements 

Initial  measurements  of  the  radiation  produced  by  the 
interaction  of  the  rotating  electron  beam  with  the  wiggler 
field  involved  inserting  a  small  horn  antenna  into  the  region 
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FIG  b  Radiated  power  spectrum  for  Bn*  =  DOUG,  =  6  28  cm. 
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FIG.  7.  Radiated  power  spectrum  for  50w  =  G.  /0  =  6.28  cm. 


immediately  downstream  of  the  interaction  region  and 
bringing  the  signal  out  through  a  length  of  Ka-band  (22-40 
GHz)  waveguide.  The  radiation  was  then  guided  through 
various  waveguide  filters  into  calibrated  attenuators  and  de¬ 
tectors.  The  receiving  horn  was  oriented  in  such  a  way  as  to 
pick  up  TM  waves  in  the  coaxial  conducting  boundary  sys¬ 
tem.  In  these  initial  experiments,  a  wiggler  field  amplitude  of 
1300  G  and  a  wiggler  period  of  6.28  cm  were  used  exclusive¬ 
ly.  Typical  signals  obtained  in  T  band  (91-170  GHz)  with 
and  without  the  wiggler  magnets  in  place  are  shown  in  Fig.  5, 
and  show  dramatically  that  high-frequency  radiation  is  only 
observed  when  the  wiggler  field  is  present.  Because  of  the 
difficulty  in  efficiently  coupling  radiation  out  of  the  system, 
we  have  only  been  able  to  estimate  the  total  radiated  power 
as  being  in  excess  of  200  kW.  Therefore,  a  reliable  measure¬ 
ment  of  the  electronic  efficiency  of  the  experiment  is  current¬ 
ly  unavailable.  These  results  have  been  reported  previous¬ 
ly.15 

In  these  measurements,  the  horn  has  been  located  in 
various  configurations,  and  the  radiated  power  has  been  ob¬ 
served  to  be  greatest  in  the  direction  of  the  electron  orbits,  as 
expected  if  the  radiation  is  scattered  in  the  forward  direction 
as  in  a  linear  FEL.  If  the  horn  is  moved  to  detect  radiation  in 
the  opposite  direction,  the  observed  radiation  is  down  by 
more  than  10  dB.  Waveguide  cutoff  filters  were  used  to  ob¬ 
tain  a  rough  idea  of  the  frequency  content  of  the  radiation. 
From  these  measurements  it  was  determined  that  most  of 
the  radiated  power  was  within  N  band  (74-140  GHz). 
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FIG.  9.  Radiated  power  spectrum  for  =  450  G,  /„  =  6.28  cm.  with  in¬ 
ner  magnets  and  inner  conducting  boundary  removed. 


Detailed  measurements  of  the  radiation  spectra  have 
been  obtained  for  a  variety  of  wiggler  configurations  using 
the  grating  spectrometer  described  previously.  These  config¬ 
urations  include  (a)  the  1300  G,  6.28  cm  period  wiggler  first 
studies  (Fig.  3,  bottom  magnet  array),  (b)  an  800  G,  6.28  cm 
wiggler  (same  array,  every  other  magnet  removed),  (c)  a  1000 
G,  3. 14  cm  wiggler  (middle  magnet  array),  and  (d)  a  430  G, 
6.28  cm  wiggler  in  which  the  inner  magnets  and  the  inner 
conducting  boundary  were  removed  entirely.  Plots  of  detec¬ 
tor  output  versus  frequency  are  shown  for  these  configura¬ 
tions  in  Figs.  6-9,  respectively. 

Several  features  of  these  spectra  are  worthy  of  mention. 
First,  the  radiation  observed  in  all  cases  is  very  narrow  band, 
approaching  the  resolving  power  of  the  spectrometer.  Sec¬ 
ond,  the  center  frequency  for  the  3. 14  cm  wiggler  (175  GHz) 
is  almost  exactly  double  that  observed  for  the  6.28  cm 
wiggler  (88  GHz),  and  both  frequencies  are  in  close  agree¬ 
ment  with  the  predictions  of  theory  summarized  in  Table  I. 
This  last  result  may  be  just  a  chance  occurrence,  as  the  pre¬ 
dictions  of  the  theory  are  quite  sensitive  to  the  values  chosen 
for  the  electron  beam  energy  and  the  applied  magnetic  field. 

Although  the  detector  used  for  these  measurements  has 
a  sensitivity  that  does  riot  vary  substantially  over  a  frequency 
range  of  10  GHz  or  so,  the  detector  sensitivity  does  fall 
steadily  as  the  frequency  of  the  detected  signal  is  increased. 
As  a  result  of  this  trend  and  the  fact  that  we  do  not  have  a 
calibration  source  at  175  GHz,  even  a  comparative  estimate 
of  the  power  radiated  at  this  frequency  is  not  possible.  If  we 
extrapolate  our  calibration  of  the  detector  response  versus 
frequency  from  our  data  in  the  range  26-1 35  GHz,  where  we 
do  have  calibration  sources,  then  the  radiated  power  at  175 
GHz  appears  to  be  roughly  comparable  to  that  at  88  GHz. 

IV.  CONCLUSIONS 

The  experimental  studies  indicate  that  the  production 
of  millimeter  wave  radiation  by  the  interaction  of  a  rotating 
electron  beam  with  an  azimuthally  periodic  wiggler  field  has 
been  achieved  in  a  manner  consistent  with  theoretical  expec¬ 
tations.  The  agreement  between  the  measured  radiation  fre¬ 
quencies  and  predicted  values  given  in  Table  I  is  excellent, 
but  further  work  will  be  need  in  order  to  determine  if  radi¬ 
ation  is  also  produced  at  the  other  frequencies  indicated. 
Although  the  efficiency  of  con  version  of  electron  beam  ener¬ 
gy  to  radiation  is  currently  low  (less  than  1  %),  it  is  not  at  all 
clear  how  efficiently  the  radiation  is  currently  being  coupled 
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out  of  the  system.  In  addition,  no  attempts  have  been  made 
to  minimize  wall  losses  or  to  optimize  the  experimental  con¬ 
figuration  for  maximum  radiated  power.  For  example,  it 
may  be  possible  to  design  a  system  that  will  be  unstable  at 
only  one  frequency,  in  contrast  to  the  several  unstable  fre¬ 
quencies  predicted  for  the  present  configuration. 

We  note  that  in  order  to  achieve  frequencies  similar  to 
those  reported  above,  the  gyrotron”  would  require  magnetic 
fields  more  that  an  order  of  magnitude  higher.  Comparing 
l u,  and  rj  of  Table  I  with  those  of  a  gyrotron  operating  with 
similar  parameters,  one  finds  that  the  gyrotron  has  some¬ 
what  lower  growth  rates  co, ,  but  higher  efficiencies  i). 

Future  studies  of  this  novel  source  of  coherent  radiation 
will  include  attempts  to  generate  submillimeter  radiation  by 
reducing  the  wiggler  period  and  experiments  using  lower 
energy  and  lower  current  electron  beams  with  longer  pulse 
durations. 
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Abstract 

The  macroscopic  fluid  equations  are  used  to  study  the  effect  of  beam 
thickness  on  the  linear  growth  rates  of  a  large  orbit  gyrotron.  The  chosen 
equilibrium  models  a  rotating  annular  electron  beam  produced  by  a  magnetic 
cusp.  Two  methods  are  used  to  calculate  the  growth  rates.  The  analytic 
approach  solves  the  problem  to  first  order  in  the  ratio  of  beam  thickness 
to  average  beam  radius  and  recovers  the  zero  thickness  limit.  The  numeric 
approach  uses  complete  orthonorraal  systems  comprised  of  the  homogeneous 
wave  equation  solutions  to  calculate  the  growth  rates  for  beams  with 
arbitrary  thicknesses.  It  is  concluded  that  thick,  tenuous  beams  in 
resonant  systems  must  be  operated  near  the  cutoff  of  the  electromagnetic 
wave.  The  decrease  in  growth  rate  as  the  wavelength  decreases  is 
attributed  to  the  loss  of  synchronism  in  systems  where  the  beam  has  a 
spread  in  its  canonical  angular  momentum. 


I.  INTRODUCTION 


Microwave  generation  at  cyclotron  harmonics  from  large  orbit  gyrotrons 
has  been  studied  extensively  in  recent  years.  The  reduced  magnetic  field 
requirements  inherent  in  high  harmonic  operation  make  such  devices  good 
candidates  for  high  power,  high  frequency  applications.  Considerable  work 
has  been  done  on  the  linear  growth  rate  problem  for  the  large  orbit 
gyrotron  configuration.  Growth  rates  have  been  calculated  from  single 
particle,  °  fluid,  »  and  Vlasov  formalisms  in  a  variety  of  waveguide 
geometries.  Most  of  these  analyses  assume  that  the  rotating  electron 
layers  have  zero  radial  thickness.  In  this  paper,  the  macroscopic  fluid 
equations  are  used  to  study  the  effect  of  beam  thickness  on  the  linear 
growth  rates. 

1 2 

A  recent  experiment  has  suggested  that  an  efficient  way  to  produce  a 
moderate  energy  rotating  beam  is  to  pass  a  linearly  streaming  beam  from  a 
thermionic  diode  through  a  magnetic  cusp.^  (Other  viable  injection 
schemes  include  gyro-resonant  rf  acceleration^  and  tilted  Pierce-type 
guns.^)  The  equilibrium  selected  for  this  study  models  the  cusp-injected 
configuration.  Although  the  two  previous  fluid  models  cited  (Refs.  7  and 
8)  allow  for  finite  thickness  beams,  they  never  consider  the  case  where  the 
beam  thickness,  axial  velocity,  and  axial  electromagnetic  wavelength  are 
simultaneously  finite.  This  situation  is  inherent  in  thermionic  cusp- 
injected  systems  and  leads  to  one  of  the  main  results  of  this  paper.  In 
Section  II,  the  fluid  equilibrium  is  presented  and  the  perturbed  density 
and  velocity  terms  are  calculated.  An  analytic  solution  is  described  in 


Section  III.  The  problem  is  solved  to  first  order  in  the  ratio  of  the  beam 
thickness  to  the  average  beam  radius.  The  result  is  shown  to  be  in 
agreement  with  the  single  particle  result  in  the  appropriate  limit.  A 
numeric  solution,  valid  for  arbitrary  beam  thickness,  is  the  subject  of 
Section  IV.  The  results  of  this  study  are  summarized  in  Section  V. 


II.  THE  FLUID  MODEL 


A.  The  Equilibrium  Distribution 

The  beam  equilibrium  is  shown  in  Fig.  1.  It  is  assumed  that  the 

equilibrium  results  from  the  passage  of  a  thick,  linearly  streaming  beam 

through  a  perfect  magnetic  cusp  (of  magnitude  Bq).  The  beam  originates  on 

a  cathode  that  has  a  thickness  of  2Ar  and  an  average  radius  of  rQ.  The 

electron  beam  is  assumed  to  be  monoenergetic  and  the  self-fields  are 

neglected.  Thus,  after  the  cusp,  all  electrons  move  on  helical  orbits 

about  the  same  gyration  axis.  The  value  of  each  electron's  post-cusp  axial 

velocity  is  a  function  of  the  injection  radius.  Any  electrons  above  a 

certain  critical  radius  will  have  insufficient  energy  and  will  be 

reflected.  That  maximum  radius  is  given  by  r  =  (y  -  1)  (m  c/eB  ) 

&  J  max  o  o  o 

2 

where  YQm0c  is  the  equilibrium  energy  (see  Ref.  13). 

For  this  idealized  situation,  the  equilibrium  electron  velocities  are 
only  a  function  of  r: 


v  (r)  =  0 
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where  Qq  =  (eBo)/(moY0)  and  yq  are  Independent  of  r.  Because  vr  «  0,  the 

o 

beam  will  be  in  equilibrium  tor  any  density  profile  which  is  only  a 


function  of  r.  In  this  paper,  two  different  density  functions  are 
considered.  The  first  is  a  piecewise  cubic  model  and  is  given  by 
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where  z  =  (r  -  r  )/Ar.  The  parameter  p  is  adjustable:  0  <  p  <  1.  A 


typical  density  profile  and  density  profile  derivative  are  shown  in 


Fig.  2.  The  parameter  n  is  a  measure  of  the  volume  density  of 


particles.  The  surface  density  of  particles  is  given  by 


r  +Ar 
o 

n  -  /  n  (r)dr  . 
s  .  o 

r  -Ar 
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For  the  piecewise  cubic  model,  the  relationship  between  the  surface  and 
volume  densities  is  given  by 


ng  =  nQ( I  +  p  )Ar  . 


(3) 


When  p  =  1,  the  distribution  function  becomes  the  usual  constant  density 
profile : 


ns  =  n0(2Ar)  *  (5) 

The  piecewise  cubic  model  is  used  in  both  Che  analytic  and  Che  numeric 
problems . 

The  second  density  function,  used  only  in  the  numeric  solution,  is 
given  by: 


n  (r) 
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n  exp[-  z  /(I 
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else 


where  z  can  be  adjusted  to  achieve  a  desired  profile 

be  shown  that  nQ(r)  is  infinitely  differentiable  and 

are  continuous  everywhere.  The  extrema  of  n'(r)  can 

o 

to  be 
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Notice  that  z  ♦  -k  1  as  z  ♦  0*  In  fact,  this  distribution  approaches  the 

A 

constant  distribution  as  z  -*•  0.  It  is  more  practical,  however,  to  use  the 
first  distribution  when  approximating  the  uniform  case. 

For  the  second  distribution,  the  surface  and  volume  densities  are 
related  by 


n  =  n  Arf  z2e_z  r( l/2)U(3/2,2,z2)]  (8) 

SO 

where  U(v,p,x)  is  a  degenerate  hypergeometric  function,^  and  T(x)  is  the 

“2  -1  *  2 

Gamma  Function.  It  can  be  shown  that  z  U(3/2,2,z  )  -*■  r(3/2)  as  z  +0 

A  *2 

which  means  that  n  ♦  n  (2Ar)  as  z  >0.  Profiles  of  n  (r)  and  n'(r)  are 
so  o  o 

AAA 

shown  in  Fig.  3  for  z  =  0.125,  z  =  1.0,  and  z  =  2.0.  The  curves  are 
normalized  so  that  nQ(z)dz  =  2.  The  curves  support  the  claim  that  the 

A 

constant  density  profile  is  approached  as  z  +  0. 

The  natural  definition  of  the  beam  thickness  for  the  numerical  routine 
is  z  =  1,  because  the  density  distributions  are  identically  zero 
when  |z|  >  1.  However,  as  suggested  by  Fig.  3,  this  definition  may  not  be 
appropriate  when  comparing  results  for  different  parameters.  Three 

A 

appropriate  measures  of  the  effective  beam  thickness  are:  (1)  n  /(2Arn  ), 

s  o 

(2)  the  location  of  the  (negative)  peak  of  n^,  and  (3)  the  Full  Width  Half 
Maximum  (FWHM).  For  the  piecewise  cubic  model,  these  measures  have  the 
common  value  z  =  (1  +  p)/2.  For  the  smooth  model,  all  three  measures  are 
different.  The  first  is  given  by  Eq .  (8),  the  second  by  Eq.  (7),  and  the 
third  by 
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can  also  be  linearized  and  Fourier  analyzed  to  yield 
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In  the  derivation  of  the  above  results,  all  variables  were  written  as 


X(x,t)  =  XQ(r)  +  X^r)  exp[i(kzz  +  £8  -  tut)] 

and  |Xj(r)|  <<  |x  (r)|  was  assumed.  Equations  (10)  and  (12)  express  the 
perturbed  velocity  and  density  in  terms  of  the  equilibrium  quantities  and 
the  electromagnetic  (EM)  fields.  These  equations  are  the  starting  points 


for  both  the  analytic  and  numeric  solutions.  The  differences  in  the  two 
solutions  arise  from  the  procedures  used  to  close  these  equations. 
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III.  THE  ANALYTIC  SOLUTION 

A.  Th,1  Thin  Beam  Approximation 

The  assumption  that  yields  the  analytic  solution  is  that  the  beam 
thickness  is  small  compared  to  the  average  beam  radius.  That  is,  the  thin 
beam  must  satisfy: 

x0  =  ~r- «  1  •  03) 

o 

The  fluid  equations  are  expanded  about  the  equilibrium  radius  in  the 

coordinate  x,  defined  by  the  transformation  x  =  r/rQ  -  1.  Throughout  the 

course  of  this  analysis,  many  quantities  will  be  linearized  in  x,  resulting 

in  additional  constraints  on  the  size  of  xq.  Unfortunately,  it  turns  out 

that  the  solutions  will  be  valid  for  reasonable  values  of  x  only  In  a  few 

o 

special  cases.  Nonetheless,  this  method  yields  a  good  deal  of  information 
about  the  linear  growth  rate  problem  as  will  be  demonstrated  in  the 
remainder  of  the  section. 

Two  critical  limitations  in  the  expansion  parameter  come  from  the 

linearization  of  vz  (r).  To  first  order  in  x, 

o 

Vz  =  v2  U  ~  x(v0  /vz  )2j  =  vz  U  -  xt)  (14) 

o  o  o  o  o 

A  A 

where  v  and  v  are  the  azimuthal  and  axial  velocities  at  the  center 
6  z 

°  o  ~  r  2 

radius  r  ,  respectively,  and  t  =  (v,  /v  )  .  The  worst  case  for  this 
o  o  z 

°  °  2 

expansion  is  at  x|  =  x  ,  resulting  in  the  criteria  x  <<  (v  /v ,  )  . 

o  o  z  0 

o  o 
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Because  large  orbit  gyrotrons  are  normally  designed  to  extract  azimuthal 

A  A 

energy,  it  is  desirable  to  have  /v  as  small  as  possible,  causing  a 

o  o 

further  restriction  on  the  allowed  size  of  xQ. 

The  zero  thickness  linear  growth  rate  solution  that  we  compare  our 
results  to  assumes  a  resonance  condition,  for  which  the  frequencies  of  the 
beam  mode  and  an  empty  waveguide  mode  coincide.  It  also  requires  a 
synchronous  condition:  the  frequency  of  the  actual  electromagnetic  mode  is 
assumed  to  be  near  the  frequency  of  the  beam  mode.  The  analysis  in  this 
section  deals  only  with  the  synchronous  case  and  the  empty  waveguide  mode 
does  not  enter  into  the  picture  directly.  Still,  the  resonant  limit  will 
be  recovered  for  the  proper  axial  wave  numbers  (kz). 

The  concept  of  synchronism  is  somewhat  ambiguous  for  the  thin  beam 
model,  because  the  phase  velocity  of  the  beam  perturbation  is  in  general  a 
function  of  radius: 


w  (r)  =  iil  +  k  v  (r)  .  (15) 

o  o  z  z 

o 


For  simplicity  in  the  analysis,  it  Is  assumed  that  the  synchronous  layer 
corresponds  to  r  =  rQ.  This  assumption  does  not  affect  the  analysis  in  any 
way,  it  simply  serves  as  a  reference  point  when  comparing  this  method  to 
others.  For  i^(r)  =  u  -  tuo(r),  the  linear  expansion  is 


♦*(r>  =  iA  +  xkzr0^o  ^  =  +  (16) 

V 

z 

o 

A  A  A 

where  to  =  w  -  w  (r  )  and  T  =  k  r  v  /v  .  Thus,  the  synchronous  limit  is 
£  o  o  z  o  d  z 

o  o 

A 

defined  by  ^  ♦  0*  The  expansion  of  ^  requires  that 

£  i 


A  *  A  1  A  A 

«  (v„  /va  KkO  •  It  was  noted  previously  that  v  /v  is 

o  z  0  £  o  z  o  J  2  0 

oo  o  o 

designed  to  be  small.  Furthermore,  it  is  known  that  i|>  <<  SI  for  tenuous 

JC>  O 

beams. ^  Although  the  method  considered  in  this  section  is  valid  for 

A 

somewhat  higher  beam  densities,  Q  will  be  small  in  general  and  the 

available  range  for  x  will  be  diminished  further.  Taking  k  <  r~*  will 

o  z  o 

lessen  the  severity  of  this  restriction.  Indeed,  at  cutoff  there  is  no 
restriction  at  all.  When  k2  =  0,  the  beam  frequency  is  independent  of  r 
and  the  electromagnetic  wave  can  be  synchronized  with  the  entire  beam. 

In  the  following  analysis,  other  expansions  in  x  are  performed. 
However,  the  linear  coefficients  of  those  expansions  are  not  inherently 
large  and  the  restrictions  imposed  by  them  are  not  calculated.  Still,  it 
will  be  demonstrated  later  that  the  solution  is  very  sensitive  to  small 
changes  in  certain  terms.  Thus,  the  error  introduced  by  one  of  these 
"accurate"  expansions  can  change  the  answer  dramatically.  The 
uncertainties  discussed  above  make  it  difficult  to  predict  the  maximum  beam 
thickness  that  yields  reliable  growth  rates.  Nonetheless,  the  method  is 
accurate  in  the  zero  thickness  limit  and  should  be  able  to  predict  the 
initial  growth  rate  trends  as  xQ  increases  from  zero. 

B.  Method  of  Solution 

The  solution  is  achieved  by  integrating  several  components  of 
Maxwell's  equations  from  rQ  -  Ar  to  rQ  +  Ar.  The  resulting  equations 
express  the  "jump"  in  the  EM  fields  across  the  beam  in  terms  of  integrals 
of  the  source  terms  and  the  EM  fields  in  the  beam  region.  In  order  to 
perform  the  integrals,  the  EM  fields  must  be  evaluated  to  first  order 
in  xQ  in  the  region  |x|  <  xQ. 


This  is  done  in  Section  1.  The  source 


integrals  are  evaluated  in  Section  2.  In  Section  3,  the  formal  solution  is 


given.  Properties  of  the  solution  are  discussed  in  Section  II. C. 


1.  The  Field  Expansions 

Outside  the  beam  region,  the  homogeneous  wave  equations  are  valid  and 
the  axial  fields  can  be  written 


E 

z 


( Vy)  — 

Vo  Ez 

Vy’V> 


z 


and 


0  <  r  <  r 


(17a) 


Tj  <  r  <  r 
+  w 


r 


cB  =  / 
z 


OoO  ~  - 

- ttt - r  (3  cB  ) 

y_Jp(y_)  r  Z 


r+Qt(y,yw)  ~  + 

*  (3  cB  )+ 


n  /  .  yu  j 

WW  r  2 


0  <  r  <  r 


r,  <  r  <  r 
+  w 


(17b) 


for 

r*  =  rQ  ±  Ar,  E*  =  E^),  (3^)*  =>  C3rc\)|r  »  ^  =  Kw/c)2  -  k2], 

y  =  r£,  y  =  y(r  ),  and  where  P  —  S  are  cross  products  of  Bessel 

^  lx, 

functions : 
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PA(a,b>  =  Jjl(a)Yjl(b) 

-  Y£(a)J£(b) 

(18a) 

Q^Ca.b)  =  J4<a)Y'(b) 

-  Va)Ji(b) 

(18b) 

R4(a,b)  =  J^ajY^b) 

-  Y;(a)Jt(b) 

(18c) 

S4(a,b)  =  J'(a)Y£(b) 

-  Y^(a)J^(b)  , 

( 18d) 

When  |x|  >  xq,  5  =  0,  so  and  can  be  found  from  Maxwell's  Equations 
given  Eq.  (17). 

In  the  region  |x|  <  xQ,  the  fields  are  expanded  about  the  values  of 
the  fields  at  the  edges  of  the  beam  (r  ).  The  piecewise  cubic  beam  density 

T 

profile  is  used  so  that  the  values  of  the  fields  and  their  derivatives  are 
continuous  at  the  beam  edges.  Thus,  the  limiting  values  of  the  homogeneous 
fields  can  be  used  in  the  expansions.  Let  x  represent  any  EM  field 
component.  For  the  fields  to  be  accurate  to  first  order  in  x  ,  it  is 

A  O 

necessary  to  approximate  k  by  a  cubic: 

xc(t)  -  lx°  -  ^  <^>4]  +  Xs  -  09) 


+  *2l^  c^>4]  -  <&>°  -  4-  x4] 


~2 


o 

where  x  *  (x  ±  X  )/2,  x  =  x(r  )*  etc.  By  considering  Taylor  series 

± 

2 

expansions  about  r  ,  it  can  be  seen  that  I X  “  Xc I  =  0(xo)  when 
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|z|  <  1.  Furthermore,  Eq.  (19)  insures  that  d\/dr  is  continuous  on  the 


same  interval. 


2.  The  Source  Integrals 

To  lowest  order  in  xq,  the  jump  in  the  perturbed  density  across  the 


beam  is  (see  Eq.  12) 


r  +Ar 
o 


TO  VAr 

O  ,e  „  t 


—  /  (rn  )dr  =  — -  ( — )  /  n  (1  +  x)v  dr 

e  '  1  ,„2e  i.o  r. 


'o  r  -Ar 
o 


V  o  r  -Ar 
l  o 


#  .  ro+Ar  to 

+  IU  <^)  /  n00  -  x  )v  dr 

°  vAr  1 


k  r  ro+Ar  TO 

+  4-^  (— )  /  n  [  1  +  x(  1  -  -4)  ]v  dr 

V  o  r  -Ar  °  *  Z1 

l  o  i 


(20a) 


Because  3.  =  -  en  v.  -  en.v  , 
1  o  1  1  o 


ro+Ar  ro+Ar  ro+Ar 

-1-  /  3  dr  -  -  1  (— )  /  n  v  dr  -  —  (— )  /  (rn  )dr 

Eoc  r  -Ar  61  0  Eo  r  -Ar  0  91  E  Eo  r  -Ar  1 

o  o  o 


(20b) 
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r  +Ar 
o 


z  TO 


r  +Ar 
o 


h  l  ..(t\  +  dr 


•o"  r  -Ar  *1 
o 


*i  \  °  VAr 


V-  f  ro+Ar  TO 

-  —  *-(-£->  /  "ot1  -  +  x)]v  dr 

*£  °  ro“Ar  h  1 


v  ,  r  +Ar 

z  k  r  .  o  Til 

-  -4-2.  (^2  )  /  nQ[ 1  ~  *(-^2  +  t  -  l)]v  di 

%  o  ro-Ar  V,  1 


r  .  ro+Ar 

—  (— 5  /  no(1  +  x)vz  dr 

o  r  -Ar  1 

o 


(20c) 


again  to  lowest  order  in  xq.  Therefore,  evaluation  of  the  source  Integrals 

ro+Ar 

requires  the  evaluation  of  integrals  of  the  form  f  n  (1  +  Fx)v.dr  where 

r  "Ar  ° 

F  Is  some  constant. 

From  the  above  comment  and  Eq.  (10),  It  can  be  seen  that  all  the 

source  integrals  can  be  expressed  in  terms  of  integrals  of  the  form 
ro+Ar 

f  n  (1  +  fx)X  dr  where  x  again  denotes  a  field  component.  The  cubic 

r  -Ar 
o 

approximation  for  x  can  be  used  to  show  that 


r0+Ar 


/  n  ( r ) ( 1  +  Fx)x(r)dr 


r  -Ar 
o 


•SU  -  y  t1  +  (!"  p)  (3  +  —)]r„  (^-)A 


o  dr' 


=2  Ll  +  4-  (1  -  p ) C 3 1  +  19p  -  5p2  -  3p3)]fxA} 


1 

to  first  order  in  xq.  Even  though  the  continuity  of  dnQ/dr  was  required  in 
the  derivation  of  Eq.  (19),  the  source  integrals  are  continuous  functions 
of  p  [their  dependence  is  given  by  Eq.  (21)].  In  the  remainder  of  this 
paper,  p  is  set  to  one  for  convenience  and  only  the  constant  distribution 
of  Eq.  (4)  is  considered. 

In  Appendix  A,  Eqs.  (10)  and  (21)  are  used  to  find  standard  integrals 
of  the  perturbed  fluid  velocity  in  terms  of  the  equilibrium  beam  quantities 
and  the  fields  (and  their  derivatives)  at  the  beam  edges.  Thus,  from 
Eqs.  (20)  and  (A.l),  the  integrals  of  the  perturbed  sources  across  the  beam 
are  known  in  terms  of  the  electromagnetic  fields  at  the  beam  edges. 


3.  The  Jump  Conditions  and  The  Formal  Solution 

Six  of  the  eight  components  of  Maxwell's  equations  contain  derivatives 
of  the  electromagnetic  fields  with  respect  to  r  and  can  be  integrated 


across  the  beam.  The  results  to  first  order  in  xQ  are: 


ro+Ar 

-iEA  +  x„(-i?°  +  iE°  +  k.rJ5°)  -  —  ^  /  (rn^dr  =  0 


r  o  r  6  z  o  z  2r  e 

o  o  r  -Ar 
o 


.  a,r 

-icB  +  x  (-icB°+  — -  E° 
0  o  0  c  z 


-icBA 


air 


x  ( — -  E°  +  k  r  cB°) 
ov  c  0  z  o  r 


ro+Ar 

1  i  r  -  , 

TFT  /  A  J8.dr 
o  r  -Ar  1 
o 


=  0 


(22a) 


ro+Af 

ic5r>  +  A  /  dr  -  0  (22b) 

o  o  r  -Ar  1 
o 


(22c) 


cSA  +  x  lc&°  -  £(-ic&°)  -  k  r  (-icS0)]  =  0 
r  oLr  6  z  o  zJ 


(2  2d) 


air 


^  +  +  *<-  1E°>  +  -r  -  o 


(22e) 


air 


+  x  [k  r  (— iE°)  -  —  ( -icB° ) ]  =  0  . 

Z  0^20  X  C  Q  J 


(22f) 


Out  of  the  six  jump  conditions,  three  contain  source  terms  and  three  do 
not.  The  three  sourceless  equations  are  linearly  dependent  as  are  the 
three  equations  with  sources.  It  can  be  shown  that 

air 

-  ( - )  x  Eq.  (22a)  +  (k  r  )  x  Eq.  (22b)  -  l  x  Eq.  (22c)  =  0 

c  z  o 

air 

(__J1)  x  Eq.  (2 2d)  +  (k,rQ)  x  Eq.  (22e)  -  l  x  Eq.  (22f)  =  0  . 

The  four  equations  (22b),  (22c),  (22e),  and  ( 2 2 f )  form  a  linearly 
independent  set  even  if  £  =  0  =  k2 .  Consequently,  they  are  the  ones  used 
in  the  analysis. 

Because  E.  and  cB.  are  given  in  terms  of  E  ,  cB  ,  and  their 

z  z 

derivatives  when  |x|  >  xQ,  Eq.  (17)  reveals  that  there  are  only  four 
unknowns.  Consequently,  there  are  the  same  number  of  equations  as  unknowns 
and  the  system  of  equations  is  well  defined.  To  simplify  the  xQ  =  0  limit, 


the  four  unknowns  are  chosen  to  be: 
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uL  =  (Sz/r)°  (23a) 
u2  =  (§2/r)A  (23b) 
u3  =  -  i(3rcSz)°  (23c) 
u4  =  -  i(3rcSz)i  .  (23d) 


Defining  the  normalized  wave  admittances 

“e  =  (r3rV±/(V*  (24a) 

and 

a*  =  (“icB  )*/(-ira  cB  )*  (24b) 

d  z  r  z 


and  using  Eqs.  (17)  and  (24)  it  follows  that 


°«  -  ty+V'Vyu)l/Wyw) 
“b  "  9»(y+’yw)/ly+VVyv)> 


ae  = 

(25) 

a~=  J£(y_)/[y_J'(y_)l  . 


The  thin  beam  geometry  parameters  are  defined  to  be  (in  the  usual 


notation) : 


o  =  (a  +  a  ) 
e  2  e  e 


ab  =  7  (ab  +  ab) 


A  1  ,  +  -v 

o  =■  -x-  (a  -  a  ) 
e  2  e  e 


A  1  ,  +  -v 

ab  =  T  (ab  "  V 


These  expressions  can  be  expanded  in  xQ  to  reveal  their  limiting  behavior 


,  y  R  yJ'AyJ  J*  (y..)  y  *  y_J'(y_) 


y.JtyJ 


y  S  VVV 
o 


°  *yoJ^yo>S 


°  yS  yoJZ(y o> 


A  r  J£(yw)  .  ,  j,  l  t  .  2  2W(,  0*2  ,  W  ,2,, 

^yoj;(y0)S  yoS  WV 


for  P  =  P4(yo.yw).  etc. 

Maxwell's  equations  can  be  used  to  find  all  the  fields  to  first  order 
in  x„  in  terms  of  u,  —  uA .  For  example,  the  average  fields  (x°)  are 
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.(lfur  kzu3>^2  (28a) 

=  ("  £kzui  +  ?  u3)/?2  (28b) 

S  =  roul  +  Xorou2  (28c> 

«  [kz(“eu!  +  “eU2>  "  1  ?  V3  +  “£\)]/52  (28d) 

_ic^0  =  (aeUl  +  °leU2)  "  ikz(abU3  +  “bVl^2  (28e) 

~LcK  =  ro(ab  +  Vb)u3  +  ro(ab  +  Xoab)uA  (28f) 

The  difference  fields  (yA)  can  be  found  from  the  average  fields  in  Eq.  (28) 

by  interchanging  u^  and  U2  and  by  interchanging  u-j  and  u^ 

(e.g.  =  r  u0  +  x  r  u, ) .  The  field  derivatives  (dy/dr)A  can  also  be 

&zo2ool  A 

found  in  terms  of  u^ — u^  directly  from  Maxwell's  Equations. 

At  this  point  the  solution  Is  complete,  however,  the  expressions  are 
too  cumbersome  to  assemble  into  one  equation.  Instead,  the  equations  are 
built-up  as  follows.  The  geometry  parameters  are  given  in  terms  of  the 
system  configuration,  the  wave  numbers  i  and  kz,  and  the  unknown 
frequency  w  [Eq.  (27)].  The  EM  fields  are  given  in  terms  of  Uj —  u^.  via 
the  geometry  parameters  [Eq.  (28)].  The  sources  integrals  can  be  found  in 
terms  of  —  u^  via  the  EM  fields  [Eqs.  (20)  and  ( A .  1 )  ]  .  Finally,  the 


jump  conditions  can  be  found  (once  again  in  terms  of  u^ —  u^)  from  the  EM 
fields  and  the  source  integrals  [Eq.  (22)]. 


C.  Properties  of  the  Solution; 

When  xq  =  0,  two  of  the  jump  conditions  yield  E^  =  0  =  u^  [Eq.  ( 22 f ) ] 
and  E^  =  0  =  u^  [Eq.  (22e)].  The  remaining  two  equations  then  can  be 

A  A 

solved  to  find  ^  and  the  ratio  u^/u^.  To  lowest  order  in  ^ ,  the 
equations  yield  the  disperson  relation 


+  —  (aV;  -  a£fl*)  =  0 
e  yq  e  0 


(29) 


With  the  aid  of  Eq.  (27),  we  can  show  that  this  expression  agrees  with  the 
single  particle  result. 

The  fluid  method  has  expanded  the  single  particle  result  in  three 
dimensions:  (1)  axial  wave  number  k_,  (2)  beam  density  n  ,  and  (3)  beam 
thickness  Ar.  Below,  the  properties  of  the  growth  rate  solutions  are 
examined  as  a  function  of  these  three  variables. 

Figure  4(a)  demonstrates  the  effect  of  beam  density  on  the  growth 

A 

rate  r  for  a  zero  thickness  beam  at  a  resonant  point  (^^  =  wr  +  ir).  The 
parameters  for  this  case  are  1  =  7,  nr  =  1  (first  radial  mode),  rQ  =  0.06  m, 

A  A 

rw  =  0.0865  m,  8g  =  0.85,  8Z  =  0.5,  and  '<z  =  0  (rw  is  adjusted  for 
°  1/3  Z°  5 

cutoff).  Because  T  «  ng  in  the  tenuous  limit,  the  plots  are  normalized 
1/3 

by  n  .  Thus,  the  "tenuous"  limit  is  valid  for  the  densities  where  the 
s 

curves  are  relatively  flat.  Consequently,  the  tenuous  limit  typically 

12-2  “ 
demands  that  n  <  10  m  in  order  for  the  lowest  term  in  i>.  to  be 
s  i 


dominant 


Self-fields  introduce  another  kind  of  tenuous  limit,  outside  of  which 
the  space-charge  effects  contribute  significantly  to  the  equilibrium  and 
(perhaps)  the  growth  rates.  Because  the  self-fields  have  been  neglected  by 
this  analysis,  there  is  an  unKnown  upper  limit  on  the  allowable  beam 
density  and  the  high  density  results  may  not  be  valid. 

Figure  4(b)  shows  T  for  the  system  parameters  of  Fig.  4(a)  midway 

between  resonant  points  (k  =  66.11  m  ^).  For  this  non-resonant 
1/2 

case,  T  «  n  in  the  tenuous  regime  and  the  plots  are  normalized 

appropriately.  The  tenuous  limit  again  appears  to  be  valid  out 
1 2  -2 

to  n  ~  10  m  . 
s 

Figure  5(a)  shows  a  plot  of  growth  rate  as  a  function  of  k  for  a  zero 

z 

thickness  beam.  The  parameters  are  i  =  11,  rQ  =  0.06  m,  rw  =  0.10  m, 

A  A  Q  A 

=  0.85,  =  0.5,  and  ng  =  10°  m~  .  The  corresponding  beam  line  and 

empty  waveguide  curves  are  plotted  in  Fig.  5(b).  The  beam  line  intercepts 
the  lowest  three  empty  waveguide  modes.  The  squares  in  Fig.  5(a)  indicate 
the  growth  rates  for  the  six  resonant  points  and  the  one  non-resonant  point 
midway  between  the  resonances  as  calculated  by  the  single  particle 
formula.  The  fluid  result  is  indicated  by  the  line.  The  growth  rates  drop 
abruptly  to  zero  at  synchronous  points  below  the  empty  waveguide  curves. 
Above  the  dispersion  curves,  the  growth  rates  drop  off  rapidly  to  small  but 
finite  values.  At  synchronous  points  sufficiently  above  the  lowest  two 
waveguide  modes,  the  growth  rates  again  become  zero,  resulting  in  five 
frequency  bands  of  instability. 

As  previously  claimed,  the  fluid  model  agrees  with  the  single  particle 
model  as  ng  -*•  0,  both  on  and  off  resonance.  It  is  interesting  to  note  that 
the  thin  beam  model's  solution  is  Independent  of  p  at  xq  = 


0.  Hence,  the 


particular  shape  of  the  density  function  is  not  important  as  xq  >  0.  This 
is  readily  seen  by  setting  xq  =  0  in  Eq.  (21)  and  recalling  that  the  source 
integrals  can  all  be  written  in  terms  of  that  equation.  The  growth  rate 
does  depend  on  the  surface  density  ng  and  the  EM  fields  at  r  ,  as  though 
the  fluid  momentum  equation  was  simply  evaluated  at  the  equilibrium  radius. 
Figure  6  shows  ^  as  a  function  of  beam  thickness  for  a  system  whose 

A 

parameters  are  t  ■  7,  n  *  1,  r  =  0.06  m,  r  =  0.09  m,  =  0.85, 

L  O  W  O 

*  8  -2  ° 

8z  =  0.5,  kz  =  0,  and  n  =  10  ra  .  Except  for  the  wall  radius,  the 

o 

parameters  of  Figs.  4(a)  and  Fig.  6  are  the  same.  The  wall  radius  is 
larger  in  the  latter  case  so  that  the  system  is  not  resonant  at  kz  =  0. 
Cutoff  is  chosen  so  that  the  entire  beam  is  synchronous  with  the  EM  wave. 
Two  curves  are  plotted  in  the  figure.  One  curve  corresponds  to  the 
solution  when  the  geometry  parameters  (a°  -  a^)  are  given  by  Eqs.  (25)  and 
(26)  (exact)  and  the  other  curve  corresponds  to  the  solution  when  the 
geometry  parameters  are  given  by  Eq.  (27)  (linear  expansion).  In  this 
example,  the  restriction  due  to  the  geometry  parameters  is  xQ  -  0.02. 

Unfortunately,  the  restriction  imposed  by  the  geometry  parameters  is 
even  more  severe  at  resonance  points.  This  can  be  seen  by  examining 
Eq.  (27).  Consider  the  solution  near  a  TE  resonant  point.  When  xQ  +  0, 
a“  +  0  as  ^  ♦  0.  Thus,  the  linear  term  becomes  the  dominant  term  for 
extremely  small  values  of  xQ.  The  solution  Is  very  sensitive  to  changes 
in  and  the  exact  and  linear  solutions  diverge  rapidly  as  xq  increases. 

In  fact,  the  thin  beam  model  is  only  useful  in  determining  the  effect 
of  beam  thickness  in  growth  rates  for  off  resonance,  cutoff  points. 

The  kz  =  0  restriction  is  due  to  the  synchronous  approxim.it ion  given  In 

A 

Eq.  (16).  For  this  special  case,  typically  varies  less  than  17  over  the 
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allowed  range  of  xq.  There  are  at  least  two  ways  to  get  more  information 
about  the  effect  of  beam  thickness  on  growth  rates.  The  first  possibility 
is  to  solve  the  equations  to  higher  order  in  xq.  The  procedure  for  doing 
this  is  straightward,  but  the  amount  of  algebraic  work  increases 
tremendously.  For  example,  in  order  to  expand  the  fields  to  second  order 
in  x  ,  they  would  have  to  be  expanded  to  fifth  order  in  z.  The  second 
order  equilibrium  density  function  would  also  be  more  complicated,  at  least 
until  p  =  1  could  be  justified.  The  second  approach  is  a  numerical  one  and 
is  the  subject  of  Section  IV. 


r‘ 


D.  Solution  of  the  Coaxial  Case 

The  only  computational  difference  between  the  cases  of  cylindrical 
waveguides  with  and  without  inner  conductors  is  in  the  calculation 
of  ag  and  a^.  The  axial  fields  between  the  inner  conductor  and  the  beam 
can  be  written 


E 
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Pjty.yp 


(30a) 


and 
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where  r^  <  r  <  r_,  and  so 
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and 
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The  new  average  and  difference  geometry  parameters  can  be  calculated  for 
the  inner  conductor  case  and  the  remainder  of  the  algorithm  does  not 
change.  The  linear  expansions  are 


.  y  R,  y  R-  P_  y  R.  y  R 

o  1  ,o  ,  o  2.  ,  ,  3  o  1  ,  o  2v 

a  =  7  (— +  — )  +  x  (  v  -x— )  ( — x —  +  — — ) 

eZP  P  irPPP  P 

*1  *2  12  1  2 


<*“-(- 


’P1P2 


>  +  *„t' 


2  2 

-  y„ 


.  i  [(i!l)2  + 


(32) 


o  1  ,  Q1  A  Q2  ,  ^  S3  w„2  2w  Q1  a  Q1  , 

«.  =  7  ( — —  +  — — )  +  x  ( —  ■  „--)U  -  y  )( — x-  +  — — > 

yoSi  yoS2  <S1S2  yoSl  yoS2 


%  -  <-A-> +  *„l>  -  i  »2  -  ^)K^-)2  - 

'*1*^  *-s.  y-s- 


o  1 


o  2 


for  =  P^(y+»yw).  P2  =  Vy-’yi>*  and  p3  =  Pt(yw’7i>*  etc* 

Figure  7  demonstrates  the  effect  of  axial  wavelength  on  ¥  for  a 
coaxial  system  with  &  *  7,  n  =  1,  rQ  =  0.075  m,  =  0.0626  m,  rw  =*  0.10  m 
d  =  0.9,  =  0.4,  and  n  =  10^  m  “. 

iJ  ’  rw  ’  P 


Figure  7(a)  shows  the  effect 
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of  kz  on  .  Between  resonance  points,  o)r  is  close  to  zero  as 
predicted.^  At  resonance,  •  T//3.  Outside  the  unstable  region, 

increases  rapidly,  indicating  that  the  solution  is  tracking  the  empty 
waveguide  mode.  Figure  7(b)  shows  the  effect  of  kz  on  r.  The  coaxial 
results  are  quite  similar  to  the  simple  empty  waveguide  results. 
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IV.  THE  NUMERIC  SOLUTION 

The  analytic  model  has  the  a  priori  restriction  that  Ar  «  r  .  Thus, 

o 

the  analytic  model  is  valid  from  Ar  =  0  to  some  small  value  which  is  at 

most  a  few  percent  of  the  average  beam  radius.  In  contrast,  the  numeric 

model  is  not  defined  for  Ar  *  0.  The  thickness  can  be  made  arbitrarily 

small,  but  the  computational  cost  increases  tremendously  with  decreasing 

thickness.  Furthermore,  there  is  no  a  priori  limit  on  the  maximum  beam 

thickness,  though  each  particular  problem  is  limited  by  the  cusp  field 

cutoff.  Thus,  the  domain  of  the  analytic  model  extends  from  the  zero 

thickness  limit  out  to  some  small  maximum  value  and  the  numeric  model's 

validity  ranges  from  some  practical  minimum  out  to  the  cusp  cutoff  layer. 

In  principle,  the  ranges  of  the  two  models  overlap  and  a  complete  growth 

rate  versus  thickness  picture  can  be  drawn.  Unfortunately,  in  most  cases 

the  thin  beam  limit  is  too  small  for  the  thick  beam  model  to  be  applied 

* 

there. 

The  numeric  approach  involves  solving  the  inhomogeneous  wave  equation 

A*  (v 

for  E  and  cB  .  In  both  simple  and  coaxial  waveguides  (among  others),  E 
Z  z  z 

can  be  written  as  a  linear  combination  of  the  homogeneous  transverse 
magnetic  (TM)  solutions  and  cBz  can  be  written  as  a  linear  combination  of 
the  homogeneous  transverse  electric  (TE)  solutions.  With  the  aid  of  these 
homogeneous  solutions,  the  inhomogeneous  wave  equation  problem  is  reduced 
to  an  exercise  in  finding  the  zeros  of  the  determinant  of  an  infinite 
matrix.  The  numerical  method  solves  various  finite  approximations  of  this 


matrix  and  compares  the  answers  to  determine  whether  or  not  the  truncated 
solutions  appear  to  be  converging.  The  apparent  limit  is  assumed  to  be  the 
correct  answer.  The  formal  solution  of  the  inhomogeneous  problem  is 
derived  in  Section  A.  Typical  results  are  discussed  in  Section  B.  In 
Section  C,  the  solution  for  the  coaxial  case  is  presented  and  results  are 
given. 


A.  The  Formal  Solution 

The  perturbed  velocity  is  still  given  by  Eq.  (10)  but  the  perturbed 
density  [Eq.  (12)]  has  to  be  found  in  terms  of  the  EM  fields  with  the  help 
of  Maxwell's  Equations: 
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The  perturbed  currents  can  then  be  found  as  before: 
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For  fixed  l  and  1  <  m  <  «,  we  define  the  functions 
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where  is  the  rath  zero  of  J„  and  a„„  is  the  ratfl  zero  of  J'  These 
im  l  im  x. 

functions  satisfy  the  orthonormal  conditions: 
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l  ui«3(r)'Wr)rdr  *  S. 
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V  1  <  m,  n  <  «,  where  6  Is  the  Kronecker  delta  function.  Furthermore,  it 

mn 

can  be  shown  that  for 
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the  expressions 
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are  valid,  i.e.  the  suras  converge  at  every  point  in  the  Interval  [0,rwJ. 

OO  00 

This  fact  is  true  because  {u^^}  and  { }  are  complete  in  the 

m=l  m=l 

appropriate  spaces  and  because  ^  and  are  assumed  to  have  the  required 
degree  of  "smoothness." 

We  define  the  unknown  vector  k  -  (tc^  < ^  *4^  by 
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^uJtm2  =  r^ru£ml  anc*  uira4  =  r^ru£m3^*  these  complete  orthonormal 

systems  (CONS),  the  numeric  problem  can  be  reformulated  as  a  problem  of 

00 

finding  allowable  values  of  aj's  and  the  corresponding  sequences  (e  ,B  |  . 

1  m  mJ  j 

To  do  this,  the  perturbed  sources  must  be  written  in  terms  of  <.  The 
sources  are  already  known  in  terms  of  )q  —  from  Eqs.  (33)  through  (35) 
and  can  be  written  symbolically  as 
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From  Maxwell's  equations  it  can  be  shown  that 
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Plugging  Eq.  (43)  into  Eq.  (42)  gives  Xj  —  Xg  terms  of  k  (notice 
that  =  <1  and  Xg  =  <3)*  Afterward,  J*  and  N*  can  be  found  from 
Eq.  (41).  The  net  result  is  that  the  perturbed  sources  are  known  in  terms 
of  <  and  can  be  symbolically  written 
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The  inhomogeneous  wave  equations  for  and  eft^  are 
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Multiplying  Eq.  (45a)  by  ruinl(r),  Eq.  (45b)  by  ru£n3(r),  and  integrating 


both  from  0  to  r  yields: 
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By  defining 

D  =  diagfC2  -  <3u/rw)2,  £2  -  (“u/rw)2, 

....  £2  “  (H,/r)2,  C2  -  (a  /r  )2,  ...]  ,  (47! 
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it  follows  from  Eqs.  (46)  and  (47)  that  the  eigenvalues  are  given  by: 

det(D  -  S)(uj)  =  0  .  (48! 

A 

If  w  is  a  solution  of  Eq.  (48),  the  fields  are  given  by  the  vector 

A  =  (E. ,  B  ,  ....  E ,  B  ,  ...)T  (49 

II  in  m 

that  satisfies  S(w)X  =  D(w)X. 

The  solution  of  the  linear  growth  rate  problem  can  in  theory  be  found 
by  locating  the  solutions  of  Eq.  (48).  Everything  that  is  needed  to  find 

A 

the  matrix  S  has  been  derived  in  this  section.  When  nQ  =  0.  S  is 
identically  zero  and  the  solutions  correspond  to  the  empty  waveguide  modes 
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as  they  must.  When  nQ  t  0,  the  solution  can  only  be  found  approximately  by 
numerical  techniques. 

Solutions  are  found  as  follows.  The  D  and  S  matrices  are  truncated 
after  a  finite  number  of  rows  and  columns.  The  integrals  in  the  S  matrix 
are  approximated  by  Simpson's  extended  rule.  The  determinant  of  the  matrix 
is  found  by  Gaussian  elimination  with  full  pivoting.  The  zeros  of  the 
determinant  are  found  by  recursive  linear  extrapolation  starting  with  the 
single  particle  results.  Although  the  thick  beam  model  is  formally  valid 
without  any  restrictions  (except  for  self-field  considerations),  it  is 
computationally  practical  only  to  look  for  solutions  near  the  tenuous, 
resonant  limit. 

B.  Properties  of  the  Solution 

Numerical  solutions  of  f  as  a  function  of  beam  thickness  are  shown  in 
Fig.  8  for  a  system  with  l  =  7,  nr  =  1  (TE) ,  rQ  =  0.06  m,  rw  =  0.0865  m, 

a  *  Q  _  O  A 

=  0.85,  $z  =  0.5,  ng  =  10  m  ,  kz  =  0,  and  z  =  2.0,  the  equilibrium 
o  o  ' 

defined  in  Eq.  (6).  For  these  parameters,  the  cusp  cutoff  limits  the  beam 
thickness  to  xQ  <  0.16.  The  solutions  are  plotted  for  several  different 
finite  approximations  of  the  D  and  S  matrices  in  order  to  demonstrate 
convergence  problems.  As  the  beam  thickness  decreases,  the  growth  rate 
solution  for  a  fixed  number  of  eigenvectors  starts  to  decrease.  The 
solution  eventually  begins  to  increase  and  ultimately  heads  toward  infinity 
as  the  beam  thickness  approaches  zero.  The  tendency  of  the  growth  rates  to 
eventually  increase  as  the  beam  thickness  decreases  Is  a  result  of  the 
truncation  of  the  D  and  S  matrices  and  Is  not  a  property  of  the  exact 
solution.  As  the  number  of  eigenvectors  used  Increases  so  does  the  region 


of  convergence.  When  50  eigenvectors  are  used,  the  result  appears  to  have 

converged  all  the  way  down  to  xQ  ~  0.06.  When  xQ  >  0.13,  20  eigenvectors 

are  sufficient  to  get  an  accurate  answer.  More  eigenvectors  are  needed  for 

convergence  as  the  beam  thickness  decreases  because  the  density 

distribution  becomes  more  and  more  like  a  delta  function  distribution. 

Thus,  it  is  the  truncation  of  the  infinite  matrices  that  puts  a  practical 

limit  on  the  allowed  values  of  x  in  the  thick  beam  model.  In  the 

o 

remainder  of  the  figures  in  this  section,  the  curves  are  plotted  for  a 
fixed  number  of  eigenvectors  (<  50),  and  regions  where  the  solutions  have 
not  converged  are  indicated  by  dashed  lines. 

In  Fig.  9,  T  is  plotted  as  a  function  of  the  effective  beam  thickness 
for  three  different  beam  profiles.  The  beam  parameters  are  still  those  of 
Fig.  8;  the  other  two  curves  correspond  to  2  =  1  and  p  =  0.95.  The 
effective  thickness  for  the  three  models  as  defined  by  Eqs.  (3)  and  (8) 

A  A 

are:  0.383  x  for  2  =  2.0,  0.603  x  for  2  =  1.0,  and  0.975  x 
oo  o 

for  p  =  0.95.  The  growth  rate  is  in  reasonably  good  agreement  for  the 
various  models.  Unfortunately,  convergence  problems  make  it  impossible  to 
determine  if  the  solutions  converge  to  the  same  point  as  xQ  -*■  0  (as 
predicted  by  the  thin  beam  model).  Because  an  effective  thickness  is  used, 
the  differences  between  the  curves  are  attributed  mainly  to  the  shapes  of 
the  density  profiles.  The  curves  predict  less  than  a  1%  shift  in  r  over 
the  whole  range  of  beam  thicknesses. 

The  three  models  (single  particle,  analytic,  and  numeric)  are  compared 
in  Fig.  10  for  the  parameters  of  Fig.  8.  (except  that  the  p  =  0.95  curve  is 
shown).  The  effective  thickness  is  used  for  the  numeric  model.  The 
allowed  range  of  x  for  the  analytic  model  is  not  well  defined  because  the 


solution  is  at  a  resonant  point.  Nonetheless,  there  is  excellent  agreement 
of  the  two  models  for  r. 

The  agreement  of  the  two  fluid  models  should  not  be  taken  for 
granted.  Although  the  models  start  with  the  same  basic  equations,  the 
approaches  taken  are  quite  different.  For  example,  the  thin  beam  model  is 
essentially  an  expansion  about  the  zero  thickness  dispersion  relation 
[Eq.  (29)]  while  the  numeric  model  is  an  expansion  about  an  empty  waveguide 
solution  (the  matrix  D) .  The  two  methods  should  only  be  expected  to  agree 
perfectly  when  the  analytic  models'  Taylor  series  expansion  is  carried  out 
to  all  orders  and  the  numeric  models'  D  and  S  matrices  are  solved  exactly. 

In  Fig.  11,  the  growth  rate  is  plotted  as  a  function  of  beam  thickness 
for  two  systems  with  different  wall  radii.  The  common  parameters  are  a  =  7, 
n  =  1,  r  =  0.06  m,  fL  =  0.85,  B  =  0.5,  n  =  108  m-2,  and  z  =  2.0. 

1  O  u  Z  S 

o  o 

The  axial  wave  numbers  are  adjusted  so  that  both  solutions  correspond  to 

resonance  points.  The  first  curve  has  kz  =  0  and  rw  =  0.0865  m,  and  the 

second  curve  has  k  =  1.165  m-1  and  r  =  0.0860  m.  For  the  first  time  a 
z  w 

solution  has  been  plotted  for  a  system  having  k  ,  v  ,  and  x  all  nonzero 

z  z  o 

o 

simultaneously.  Consequently,  this  is  the  first  result  for  a  system  whose 
E  layer  is  only  approximately  in  synchronism  with  the  EM  wave  (for  the 
first  time  is  a  function  of  r).  The  finite  kz  growth  rate  is  a  rapidly 
decreasing  function  of  beam  thickness.  For  example,  the  finite  k^  growth 
rate  is  roughly  one  half  of  the  zero  thickness  growth  rate  (indicated  by 
the  square)  at  a  relative  thickness  of  7.5%. 

Because  the  two  cases  in  Fig.  11  are  very  similar  physically  (the  wall 
radius  changes  by  less  than  0.6%),  the  difference  in  the  curves  Is 
attributed  mainly  to  the  loss  of  complete  synchronism.  Tills  quasi- 


synchronous  effect  can  be  scaled  to  higher  densities  as  follows.  An 
estimate  of  the  relative  deviation  from  synchronism  at  the  outer  beam  edge 
is  [from  Eq.  (15)] 
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so  in  the  tenuous  limit- 


,  1/3. 

(k  /n  )  . 
z  s 


(50b) 


Thus,  for  every  three  orders  of  magnitude  increase  in  beam  density,  the 

growth  rate  curve  will  be  similar  for  a  one  order  magnitude  decrease  in 

axial  wavelength.  The  axial  wavelength  for  the  finite  kz  case  is  -  5.4  m. 

6  -2 

For  a  beam  density  of  n  =  10  m  ,  the  same  relative  decrease  in  r  should 

s 

10  -2 

be  found  for  A  ~  25  m.  Likewise,  for  n  *  10  ra  ,  the  appropriate  A 
z  s  z 

is  ~  1.16  m.  Normalized  growth  rates  at  isolated  points  are  indicated  in 

6  -2  10  '  —2 

Fig.  11  for  n  =  10  m  (diamonds)  and  n  *  10  m  (circles).  The  axial 
s  s 

wave  numbers  are  scaled  accordingly.  The  results  support  this  simple 

10  —2 

scaling  law.  It  should  be  noted  that  the  xQ  =  0.04  point  for  ng  =  10  m 
had  not  completely  converged  after  50  eigenvectors  and  the  actual  result  is 
(probably)  closer  to  the  curve  than  indicated. 


C.  Solution  of  the  Coaxial  Case 

The  only  difference  between  the  solutions  of  the  coaxial  case  and  the 
simple  empty  waveguide  case  is  in  the  CONS  used.  The  two  CONS  for  the 
r e  ! :< i  1 1  rase  are  given  by 
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for  fixed  l  and  1  <  m  <  ®,  and  for  y  =  r/rw  and  *  r^/r^  Now  a£m^£m^ 
is  the  mth  zero  of  (o^^  »a£m) ^£rayi ’ 1  *  E9uations  (37)  and  (38) 
are  valid  for  the  functions  defined  in  Eq.  (51)  provided  that  the  lower 
integration  limit  is  set  to  ri  [<2  and  <4  are  still  given  by  Eq.  (AO)]. 

Growth  rates  for  coaxial  cases  with  finite  axial  wave  numbers  are 
shown  in  Fig.  12.  The  parameters  are  Jt  »  7,  n  =  1,  r^  »  0.0626  m, 

rQ  -  0.07509  m,  =  o.l  m,  00  =  0.9,  3Z  -  0.A,  k ^  =  0.3  m 

8-2  0  0  - 

and  n  ==  10°  m  .  The  curves  are  drawn  for  z  =  1.0  and  z  =  2.0.  The 
s  • 

curves  are  plotted  as  a  function  of  the  relative  beam  thickness  and 
the  kz  =  0  (r  =  0.07500  m)  growth  rate  is  shown  as  a  reference.  The  two 
finite  kz  growth  rates  are  in  good  agreement  with  each  other.  As  usual. 


the  curves  drop  off  rapidly  as  the  beam  thickness  increases. 


V.  CONCLUSIONS 


I 

In  summary,  we  have  carried  out  a  linear  fluid  analysis  that  is 
applicable  to  cusp-injected  beams  with  finite  radial  thickness.  Two 
approaches  were  taken.  In  the  analytic  approach,  we  expanded  the  fluid 
equations  in  the  ratio  Ar/rQ  and  recovered  the  single  particle,  tenuous 
beam  result.  Futhermore,  we  noted  that  the  analytic  result  was  valid  off 
resonance  and  at  somewhat  larger  beam  densities  (than  the  single  particle 
result).  The  latter  extension  was  useful  in  helping  to  quantify  the  limits 
of  the  tenuous  regime.  Unfortunately,  we  discovered  that  the  analytic 
result  has  only  limited  usefulness  in  the  analysis  of  thick  beam 
configurations. 

In  the  numeric  approach,  we  solved  the  Inhomogeneous  wave  equations 
via  expansions  in  the  homogeneous  solutions.  Because  of  computational 
limitations,  results  were  only  displayed  for  resonant  systems  with  tenuous 
beams.  We  demonstrated  that  resonant  systems  with  thick,  tenuous  beams 
cannot  be  operated  far  from  cutoff  without  a  substantial  reduction  in  the 
linear  growth  rate.  We  attributed  the  reduced  growth  rates  to  the  loss  of 
synchronism  inherent  in  systems  that  have  a  spread  in  canonical  angular 
momentum  (manifested  by  a  finite  Av2  in  our  model). 

In  our  analyses  we  focused  on  two  particular  conducting  wall 
geometries:  circular  and  coaxial  waveguide  systems.  Extending  the 
analyses  to  other  waveguide  geometries  is  straightforward.  In  the  analytic 
model,  only  the  geometry  parameters  (normalized  wave  admittances)  need  to 
be  adjusted.  For  the  numeric  model,  results  can  be  obtained  if  the 
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APPENDIX  A 

Integrals  of  the  Perturbed  Fluid  Velocity 


Equations  (10)  and  (21)  are  used  many  times  to  obtain  the  following 
results : 
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Figure  Captions 


FIG.  1.  The  fluid  equilibrium. 


FIG.  2.  The  piecewise  cubic  equilibrium  density  profile:  (a)  n  (z)  and 

o 

(b)  n^(z)  for  p  =  0  (0),  0.5  (V),  and  0.95  (A).  For  each 

curve,  n  =  2Ar. 
s 

FIG.  3.  The  smooth  equilibrium  density  profile:  (a)  n  (z)  and  (b)  n'(z) 

o  o 

A 

for  z  =  0.125  (A),  1.0  (V),  and  2.0  (Q) .  For  each  curve, 

n  =  2Ar. 
s 

FIG.  4.  The  dependence  of  the  growth  rate  T  on  beam  density:  (a)  at 

resonance  (k  =  0)  and  (b)  off  resonance  (k  =  66.11  ra~*).  The 
z  z 

common  parameters  are  l  =  7,  r  *  0.06  ra,  r  =  0.0865  m, 

o  w 

A  A 

8Q  =  0.85,  8  =*  0.5,  and  x  =  0. 

0  z  o 

o  o 

FIG.  5.  Hie  dependence  of  the  growth  rate  T  on  axial  wave  number:  (a)  F 

along  a  beam  line  and  (b)  the  corresponding  dispersion  curves. 

Parameters:  l  =  11,  n  =  1-3,  n  =  10®  ra-2,  r  =  0.06  m, 

r  s  o 

A  A 

r  =  0.10  ra,  8q  =*  0.85,  8  =  0.5,  and  x  =0.  Single  particle 

W  Q  z  o 

o  .  o 

results  are  indicated  by  □  . 

FIG.  6.  The  dependence  of  the  growth  rate  r  on  beam  thickness. 

8  -2 

Parameters:  £  =  7,  n  *»  1,  n  =  10  ra  ,  r  *  0.06  m, 

r  s  *  o 

A  A 

r  =  0.09  ra,  8„  =  0.85,  8  =  0.5,  and  k  =0.  Linear 

w  9  z  '  z 

o  o 

expansion  V,  exact  solution  A. 


FIG.  7. 


The  dependence  of  t|»^  on  the  axial  wave  number  for  a  coaxial 

system:  (a)  and  (b)  T.  Parameters:  t  -  7,  n  ■  1, 

n  **  1011  ra"2,  r  =>  0.075  ra,  r,  =  0.0626,  r  =  0.10  ra,  ft  -  0.9. 

s  o  i  w  0  * 

*  o 

8  *  0.4,  and  x  = 

z  o 

o 


0 


FIG.  8 


The  convergence  of  the  truncated  solutions  for  f:  20  (A),  30  (V) 


FIG.  9. 


FIG.  10. 


FIG.  11. 


FIG.  12. 


40  (C)»  and  50  (O)  eigenvectors.  Parameters:  A  ■  7,  n  =  1, 

rQ  =  0.06  m,  ry  =  0.0865  m,  30  =  0.85,  fl  =  0.5,  n  =  108  m-2, 

k  =0,  and  z  =  2.0. 
z 

The  effect  of  the  density  profile  on  T.  The  common  parameters 
are  l  =  7,  n  =  1,  n  =  108  m~2,  r  =  0.06  m,  rIt  =  0.0865  m, 

r  &  o  w 

A  A 

0Q  =  0.85,  g  =  0.5,  and  k  =0.  Models:  □  single  particle, 
o  z  z 

0  0A 

Op-  0.95,  A  z  =  1.0,  and  V  z  =  2.0. 

Comparison  of  the  three  methods  of  solution  for  r.  The  common 
parameters  are  those  of  Fig.  9.  Models:  □  single  particle, 

O  analytic,  and  Q  numeric:  p  =  0.95. 

The  effect  of  axial  wavelength  on  the  growth  rate  r.  Common 

A  A 

parameters:  t  »  7,  n  ■  1,  rQ  -  0.06  m,  0Q  =  0.85,  8  =  0.5, 

o  zo 

*  _  O  O 

and  z  =  2.0.  Models:  □  single  particle:  A  ns  =  10  m  , 

r  ,  =  0.0865  m,  and  k  =  0;  V  n  =  108  m  2,  r  =  0.0860  m, 
w  z  s  w 

and  kz  =  1.16  m_1;<)  ng  =  106  m-2,  rw  =  0.0864  m, 

and  kz  =  0.251  ra_1;  andO  ng  =  1010  m-2,  rw  =  0.0843  m, 

and  k  ■  5.41  m 
z 

The  effect  of  axial  wavelength  on  the  growth  rate  of  a  coaxial 

8  -2 

system.  Common  parameters:  £  *  7,  nr  *  1,  n  “  10  tn  , 

rQ  =  0.075  m,  ri  =  0.063  m,  ry  =  0.10  m,  0Q  =  0.9, 

«  ° 
and  $z  =  0.4.  Models:  O  single  particle, 0  kz  =  0, 
o 

A  k  =  0.3  ra  *  and  z  =  1.0,  and  7  k  =  0.3  m  *  and  z  =  2.0. 
z  z 
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